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Abstract 


The evolution of cooperation in cellular groups is threatened by lineages of cheaters that 
proliferate at the expense of the group. These cell lineages occur within microbial communi- 
ties, and multicellular organisms in the form of tumours and cancer. In contrast to an earlier 
study, here we show how the evolution of pleiotropic genetic architectures—which link the 
expression of cooperative and private traits—can protect against cheater lineages and allow 
cooperation to evolve. We develop an age-structured model of cellular groups and show 
that cooperation breaks down more slowly within groups that tie expression to a private trait 
than in groups that do not. We then show that this results in group selection for pleiotropy, 
which strongly promotes cooperation by limiting the emergence of cheater lineages. These 
results predict that pleiotropy will rapidly evolve, so long as groups persist long enough for 
cheater lineages to threaten cooperation. Our results hold when pleiotropic links can be 
undermined by mutations, when pleiotropy is itself costly, and in mixed-genotype groups 
such as those that occur in microbes. Finally, we consider features of multicellular organ- 
isms—a germ line and delayed reproductive maturity—and show that pleiotropy is again 
predicted to be important for maintaining cooperation. The study of cancer in multicellular 
organisms provides the best evidence for pleiotropic constraints, where abberant cell prolif- 
eration is linked to apoptosis, senescence, and terminal differentiation. Alongside develop- 
ment from a single cell, we propose that the evolution of pleiotropic constraints has been 
critical for cooperation in many cellular groups. 


Introduction 


There is widespread cooperation in cellular groups where cells invest in costly traits that bene- 
fit all cells in the vicinity, such as bacteria that secrete an extracellular enzyme to digests nutri- 
ents or the more complex coordinated phenotypes of multicellular organisms. Cooperative 
traits can require that cells forego their own reproductive interests in favour of the reproduc- 
tive interests of the group as a whole [1,2]. This effect, in turn, can lead to the evolution of non- 
cooperative lineages—sometimes known as “cheaters”—that make use of collective benefits 
without investing in them and threaten cooperative function [3]. 
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The potential for cheater lineages is well documented in microbes. Mutants lacking a range 
of cooperative traits have been shown to outcompete wild-type cells [4-10] and occur in the 
field and clinic [11-13]. Such observations beg the question of how cooperation persists over 
evolutionary time. A key explanation is that many cellular groups, both in microbes and multi- 
cellular organsims, are recently derived from a single cell (clonal) [1,2,14,15]. In the terminology 
of sociobiology, this leads to high relatedness and kin selection, which is a major driver of coop- 
eration across many systems [16,17]. The argument is that, when cell groups are clonal, interac- 
tions between cooperative and cheater genotypes are prevented, which allows cooperative 
genotypes to prosper as cheater genotypes lose the shared benefits of cooperation. While kin 
selection is undoubtedly important [1,2,14,15], this explanation neglects a key feature of the 
biology of cellular groups: Even in a group founded from a single cell, cooperation can still 
break down due to the emergence of mutant noncooperators from within [18-20] (Fig 1). 
Indeed, with nonzero mutation rates, the question is when, not if, these lineages will emerge. 

Given the potential for cheater lineages, it has been suggested that genetic architecture can 
evolve to help stabilise cooperation. When the expression of a cooperative trait is linked to a 
private trait that helps a cell to survive or divide, mutations that ablate cooperation can also 
ablate the private trait and, thereby, stop the evolution of cheater lineages. This pleiotropic 
linkage of cooperative traits and private (personally benefitial) traits has been identified experi- 
mentally to be a mechanism that helps to promote cooperation in multiple microbial species 
[21-25]. In addition, an agent-based model of microbes found that links between metabolic 
and secretion genes helped to stabilise cooperation [26], and a theoretical study suggested that 
pleiotropy can promote niche construction, which is related to cooperation [27]. These studies 
raise the possibility that the evolution of pleiotropy might be a general mechanism to promote 
cooperation in cellular groups [28]. 

However, a recent theoretical paper argued broadly against the idea that pleiotropy is an 
explanation for the evolution of cooperation [29]. In particular, the authors argued that pleiot- 
ropy only evolves under conditions when kin selection is already operating to stabilise cooper- 
ation (Fig 2 in [29]), with, at best, a very minor impact on the evolved level of cooperation (i.e., 
seen in S14 and S16 Figs but not Fig 2 in [29]). As such, they concluded “Pleiotropy does not 
help stabilise cooperation over evolutionary time—cooperation is only favoured in the region 
where Hamilton’s rule is satisfied because of indirect fitness benefits.” 

As we disuss in detail in the Supporting information (S1 Text), a limitation of this model is 
that it did not explicitly capture group-level birth and death events or allow groups to develop 
for long enough to see the importance of pleiotropy for stabilising the evolution of coopera- 
tion. The time allowed for groups to develop is important because the problems with cheater 
lineages only becomes apparent as groups age (Fig 1). Pleiotropy only becomes subject to sig- 
nificant between-group selection, therefore, in longer-lived groups when cheater mutants have 
time to threaten the group. We show in the Supporting information that increasing the length 
of time that groups live for increases the levels of cooperation that evolve via pleiotropy in the 
model of [29]. However, problematic assumptions such as unbounded explonential growth 
prevented us from exploring this effect further (S1 Text). We, therefore, decided to develop a 
novel age-structured multilevel selection model for the evolution of cooperation in cellular 
groups, including both microbes and multicellular organisms. Our model predicts that pleiot- 
ropy is a powerful way to promote the evolution of cellular cooperation. 


Results 


Weare interested in understanding how multicellular groups founded by cells with pleiotropic 
constraints function as compared to groups founded by otherwise similar cells that lack these 
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Fig 1. The problem of cooperation and the evolution of pleiotropy. (A) Mutation-selection dynamics can 
undermine clonal multicellular groups. Mutation of cooperative cells (blue) can generate noncooperative cells (dark 
grey) that do not pay the costs of cooperation but reap its rewards. These cheater lineages can spread within 
multicellular groups because they divide more rapidly than wild-type cells. Invasion of spontaneous mutant lineages 
can lead to a breakdown in group function (distorted shape). (B) Pleiotropy promotes cooperation in our model 
because it leads to between-group selection on the rate at which cooperation breaks down within groups as they age. 
Beginning with groups founded by 1 pleiotropic and 2 nonpleiotropic lineages, we see that groups founded by the 
nonpleiotropy lineages loose function and are eventually replaced by the pleiotropic lineage. Cells are coloured in 3 
parts according to whether they display cooperation (blue), a privately beneficial trait (red), and pleiotropy (yellow), 
whereby the cooperative and private traits are linked. Mutations can make some traits inactive (grey). When 
cooperation is lost in a pleiotropic cell lineage (top row), the cells also lose their private trait, which stops them from 
proliferating. Group function is thereby preserved. When cooperation is lost in a nonpleiotropic lineage (second row), 
cheater lineages emerge that spread and group function is lost. As a result, pleiotropic groups thrive relative to 
nonpleiotropic ones and seed more groups, giving rise to between-group selection for pleiotropy. 


https://doi.org/10.1371/journal.pbio.3001626.g001 





constraints. We follow the effects of pleiotropic links between cooperative traits (that benefit 
the whole group) and private traits (that benefit the individual cell that carries them) on multi- 
level selection dynamics using an age-structured modelling approach (see Methods). A group 
in our model is intended to capture a group of microbes or a proto-multicellular organism, 
which lacks the division between germ and soma. Groups start from a single cell and display 
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logistic growth up to a carrying capacity K, which defines the size of the group at maturity. A 
second parameter, A, determines the expected life span of a group. This is important because it 
impacts on the amount of cell turnover that is expected after a group reaches reproductive 
maturity. Such cell turnover can be major contributor to the number of cell divisions within a 
multicellular group. For example, high rates of cell turnover occur in bacteria, which com- 
monly live attached to surfaces in structures known as biofilms where dispersing cells are 
replaced by dividing cells below them [1]. It is also common In multicellular organisms: Mem- 
bers of the genus Hydra (Fig 1) can live for several years, while their epithelial cells are esti- 
mated to turnover every few days [30]. 

The fact that groups start from a single cell in our model ensures high relatedness and 
strong kin selection, which is consistent with microbes that grow in clonal patches [1] and the 
biology of multicellular organisms [20]. However, we later reduce this within-group related- 
ness to study its effects on pleiotropy and cooperation. To study the effects of pleiotropy on 
cooperation, our modelling has to capture the stochastic effects of mutations. For this reason, 
the heart of the model is a stochastic simulation that captures populations of cells as they grow, 
and potentially mutate to other genotypes, within a group. However, as we discuss later, we 
also need to capture the evolutionary effects of cooperation at the group level, which is done 
with partial differential equations (PDEs) that allow us to capture a large (infinite) number of 
competing groups. We hope that this novel approach—stochastic simulations embedded in 
PDEs—will prove useful to understand a wide range of traits under multilevel selection 
(Methods). 


Pleiotropy slows the breakdown of cooperation within cellular groups 


Regulatory networks, and the maps from genotype to phenotype, are often complex [31]. Evo- 
lutionary models of cooperation typically overlook this complexity and instead study optimal 
trait values, an approach known as the phenotypic gambit [32]. Here, we treat genetic architec- 
ture as a trait, like any other, that can itself evolve in response to natural selection [26,29]. We 
do this with a simplified model of pleiotropy. The definition of pleiotropy can vary between 
disciplines and authors [29,33-35], and here we mean the commonly used definition: Pleiot- 
ropy is when a single locus affects 2 or more traits [24]. Specifically, our model captures how 
mutation at a given locus affects 1 cooperative and 1 private trait (Fig 2). While there are a vast 
range of possible regulatory networks that might influence any 2 traits of interest, the impacts 
of pleiotropy can be captured by a single value ¢, which is the probability that a mutation in a 
network with an active, cooperative, and private trait will give rise to a pleiotropic effect 

(Fig 2A). 

There are, therefore, 3 traits in our model: a cooperative trait, a private trait, and a pleiot- 
ropy trait, which gives rise to 8 possible genotypes (Fig 2B). The goal of our model is to explore 
which of these 8 genotypes is favoured by natural selection in the long term due to competition 
among cells within and between groups. For a given cell, each of the traits is in either an active 
or inactive state. At the heart of the model is a tension between selection for cooperation 
within and between groups, i.e., the cooperative trait decreases a cell’s relative division rate 
within a group but brings benefits to the group function as a whole. This trait might, for exam- 
ple, represent a secreted enzyme that helps nutrient acquisition in a microbial group, or the 
suppression of cell proliferation in a simple multicellular organism to ensure proper function- 
ing [18]. By contrast, the private trait simply increases a cell’s survival rate within its group. 
This trait might represent an enzyme involved in central metabolism, for example. With plei- 
otropy, mutations that ablate the cooperative trait increase the probability that the private trait 
is also lost and vice versa [29], where ¢ determines this probability (Fig 2A). A mutation matrix 
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Fig 2. Modelling pleiotropy. (A) One way to measure the strength of pleiotropy in real-world regulatory networks is 
to compute the ratio of those mutations in the network that simultaneously impact 2 terminal phenotypic traits to the 
number of genes. We call this measure ¢. In principle, this value can be calculated for any network regulating 2 traits. 
We show illustrative networks that would generate varying levels of pleiotropy across the range of values of ¢. (B) 
Genotype-phenotype map for our mathematical model of pleiotropy. To keep things simple, we model 3 traits: a 
terminal cooperative trait, z° (blue circles), a terminal private trait, z‘ (red circles), and a pleiotropic regulation trait, z 
(yellow circles). All traits in our model can be either active (coloured as above) or inactive (grey circles). This means 
there are 8 possible genotypes, each labelled g;, where i€{1,2,. . .,8}. Although pleiotropy can be active or inactive for a 
given genotype, we use a parameter $€[0,1] to tune the strength of its effect. This allows us to use our simple 
genotype-phenotype map to model biological scenarios when pleiotropy is expected to be weak as well as strong. 
Furthermore, we assume that pleiotropy only has a functional effect when both other traits are active. Thus, although 
active pleiotropy is present in genotypes g5, 24, go, and gg, we assume it only affects genotype gg, where it influences the 
probability, ¢, that a gg cell experiences a pleiotropic mutation given that a loss-of-function mutation has deactivated 
one of its other traits. 


https://doi.org/10.1371/journal.pbio.3001626.g002 





specifies the transition probabilities between all 8 genotypes in the model as a function of the 
strength of pleiotropy (see Methods). 

Under a null model in which pleiotropy does nothing, mutations affect each trait indepen- 
dently, where mutations that cause loss-of-function in a trait occur with rate y, and gain-of- 
function trait mutations occur with rate vu, where v<1. This value reflects the fact that it is typ- 
ically easier to break trait functionality than to restore or create it, and we typically take v = 
0.01 to capture the strength of this bias. We use y = 0.0001 per generation for the base muta- 
tion rate in most analyses, which describes the probability that a trait is mutated—and function 
is lost—per cell division. This value is expected to vary widely between systems and traits and 
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is intended only as an illustration. We later perform parameter sweeps of both y and v across 
several orders of magnitude. 

We begin by following the evolutionary dynamics within a group. In our first model, each 
group is founded by a single cell, which gives 8 possible group types corresponding to the 8 
cell genotypes (Fig 2B). While groups all start their life with clonal expansion of their founder, 
mutation-selection processes mean that their genotypic composition may change through 
time as they age. We can describe this process for each of the 8 group types. The dynamics for 
groups founded by genotypes 1 to 7 are shown in the Supporting information (S1-S7 Figs), 
and we focus here on genotype 8 groups (Fig 2B), hereafter referred to as “pleiotropic coopera- 
tors,” because they capture the effects of pleiotropy on cooperation (Fig 3). Groups with these 
genotypes initially grow towards their carrying capacity by clonal expansion but, depending 
on the strength of pleiotropy, have the potential to be invaded by cheater lineages that lack the 
cooperative trait but express the private trait (genotype 4; see Fig 2A). Importantly, we see that 
the extent and rate of invasion of the cheater lineage is diminished as the strength of pleiot- 
ropy; ¢; is increased. Cheater lineages make up 25% of the group by approximately day 25 in 
groups without pleiotropy, by day 40 in groups with intermediate pleiotropy, and never (not 
before 50 days) in groups with strong pleiotropy (Fig 3B). 

The resistance to invasion by cheater lineages occurs because pleiotropy reduces the fre- 
quency with which mutations give rise to a cheater phenotype. As a result, pleiotropy is able to 
increase the level of cooperation in groups (Fig 3C). The distribution of fitness effects (DFE) of 
loss-of-function mutations helps to show why pleiotropy is an effective mechanism for limiting 
cheater cell lineages within a given group (Fig 3D), something also clear from the dos Santos 
study [29]. When the strength of pleiotropy is relatively weak, mutations to the cooperative 
trait frequently give rise to mutant descendants that have a competitive advantage over the 
cooperative cells within the group. By contrast, when the strength of pleiotropy is relatively 
strong, mutations tend to have have neutral or deleterious effects on cells because a loss of 
cooperation also comes with a loss of the private trait. In the model, we assume that the effects 
of expressing the cooperative and private trait on within-group fitness are equal and opposite 
in magnitude, which is what leads to neutrality when both are lost. Some examples suggest that 
the loss of a private trait may have a stronger negative effect, such as cell death via apoptosis 
[36-38]. Such examples may lead to a negative change in within-group fitness when both traits 
are lost. We do not consider this case explicitly here, but it is only expected to strengthen the 
ability of pleiotropy to remove potential cheater lineages and thereby improve group function. 


Pleiotropy evolves to suppress cheater lineages and promote cooperation 


Our within-group model supports the established, and intuitive, argument that pleiotropic 
links between a cooperative and private trait will help to maintain cooperation [21-25]. How- 
ever, this model simply assumes that this pleiotropy exists, rather than explaining how it 
evolved. One origin of pleiotropy is a result of natural selection on traits unrelated to coopera- 
tion [25,39]. Pleiotropy is extremely common in all genotype to phenotype maps, whether or 
not cooperative traits are involved. As a result, cooperation may become pleiotropically linked 
to private traits simply through the way that regulatory networks normally evolve. A more 
intriguing alternative is that pleiotropy evolves to promote cooperation [25,31]. However, as 
discussed above, a recent theoretical treatment of this idea argued that, despite the abilitiy of 
pleiotropy to decrease the cheater load within groups, it will not generally evolve to promote 
cooperation over evolutionary time [29]. Specifically, this earlier study found that “(1) pleiot- 
ropy does not stabilise cooperation, unless the cooperative and private traits are linked via a 
genetic architecture that cannot evolve (mutational constraint); (2) if the genetic architecture is 
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Fig 3. The evolution of cooperation within groups depends on the strength of pleiotropy ¢. Within-group 
mutation selection dynamics are shown for a group founded by a cell with genotype gg, which actively expresses a 
cooperative trait, z¢ = 1, a private trait, z¢ = 1, and a pleiotropy trait, a= 1. Growth of the group as its age, y, 
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increases, is logistic, with a carrying capacity K = 200 at which point cells continue to divide, die, and turnover (see 
Methods). Dynamics are shown from left to right for 3 different strengths of pleiotropy, ¢, where pleiotropy is absent/ 
has no effect in the left-hand side column (¢ = 0). For comparison, the vertical dashed line in (A-C) shows the point at 
which noncooperative lineages shown in orange make up 25% of the group. Pleiotropy leads to this point being 
delayed (¢ = 0.5) or prevented (¢ = 1). (A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative 
frequencies, x“(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, Z‘(y). (D) 
Distribution of fitness effects: Shown is the effect on within-group fitness of the different types of loss-of-function 
mutations that occur, which are coloured by their phenotypic effect, where blue is loss of cooperative trait, red is loss of 
the private trait, yellow is loss of pleiotropy, and brown is loss of both cooperation and private trait (due to pleiotropy). 
When pleiotropy is weak or absent, loss-of-function mutations to cooperative traits increase the fitness of cells within 
the group, and loss-of-function mutations to private traits decrease the fitness of cells within the group. When the 
strength of pleiotropy is increased, mutations to either trait tend to have pleiotropic effects, which cancel one another 
out, meaning mutant cell lineages no longer gain an advantage within the group. Formally, the fitness effect is 
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DFE = In(b¢,,,/b¢,) — In(d¢,,,/d‘,,), where b°, and b‘, are the within-group birth rates of the mutant descendant and 
ancestral wild type, respectively, and d‘,, and d‘, are the death rates of the mutant descendant and ancestral wild type, 
respectively. Parameters: s = s* = 0.95; K = 200; u = 0.0001; v = 0.01. The code required to generate this figure can be 


found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 
https://doi.org/10.1371/journal.pbio.3001626.g003 





constrained in this way, then pleiotropy favours any type of trait and not especially cooperation; 
(3) if the genetic architecture can evolve, then pleiotropy does not favour cooperation; and (4) 
there are several alternative explanations for why traits may be linked, and causality can even be 
predicted in the opposite direction, with cooperation favouring pleiotropy.” 

The authors did find that pleiotropy will often help reduce the prevalence of cheaters within 
certain groups, in line with experimental evidence in microbes [21-23]. However, they also 
found that between-group selection was not strong enough to notably increase cooperation 
across the whole population over evolutionary time (Figs $14 and S16 in [29] show a small 
increase in evolved cooperation with pleiotropy). We discuss this study in detail in the Sup- 
porting information (S1 Text), where we conclude that a different type of model will benefit 
the study of the evolution of cooperation and pleiotropy. 

We developed our model, therefore, to study whether pleiotropy will evolve as a mechanism 
to promote cooperation within multicellular groups over evolutionary time. To do this, we 
extend the model to capture how groups perform, and compete, across a wider population. In 
each group, within-group evolution occurs as just discussed (Fig 3), which, in turn, affects a 
group phenotypic trait, which we call the group’s “function.” For example, this might represent 
the ability of a bacterial strain to produce a protective biofilm, or the ability of a multicellular 
organism to coordinate its development. The key is that the group function is assumed to 
breakdown when within-group selection dynamics take hold. 

We capture group function at age y by z;(y) = Z‘,(y)z;(y), where Z;(y) is the average 
amount of cooperation in the group, and Z‘(y) is the average expression of the private trait 
(see Methods for more details). Thus, we assume that cells expressing both cooperative and 
private traits contribute fully to group function, whereas cells with either trait missing do not. 
Specifically, cells lacking only cooperation behave as cheater lineages, as just discussed, while 
cells lacking the private trait function poorly. As a result, groups founded by these cells (geno- 
types 1 to 6; Fig 2B) have a greater probability of extinction than higher functioning groups 
(genotypes 7 and 8; Fig 2B). But within-group evolutionary dynamics also means that group 
function can degrade in groups founded by cooperative cells (Fig 3). 

We study the evolution of these processes with a system of PDEs, which captures an infi- 
nitely large age-structured population of groups. We start with a population dominated by 
individuals lacking all traits (genotype 1) and allow them to evolve by numerically integrating 
within-group and between-group dynamics over time until a stable age distribution is reached. 
We then ask if pleiotropy evolves as a function of 3 key parameters: the strength of pleiotropy, 
¢, expected group life span, A, and group size at maturity, K. 

A common approach to model multilevel selection in evolutionary biology is the haystack 
model [40], which is the approach used in [29]. In the haystack model, group selection is mod- 
elled implicitly as an emergent process resulting from fitness differences between individuals 
within groups; there is no explicit consideration of group-level survival and reproduction events. 
Instead, all group phenotypes are described in terms of the set of individual phenotypes. In con- 
trast, the PDE approach we use models group-level phenotypes and considers their direct impact 
on group-level survival and reproduction events. This approach allows us to more explicitly cap- 
ture microbial groups and multicellular organisms, which have group-level traits that are impor- 
tant for group fitness but are threatened by within group competition. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001626 June 3, 2022 8/34 


PLOS BIOLOGY Pleiotropic constraints promote the evolution of cooperation in cellular groups 













Trait 
Cooperative Private Pleiotropy 
| | 2 

TT 

~~ 

S 

Ss 
: pe 
¢ 9 
f x § 
@ | 
S 2 
£ © 
(0) 

i 

x 

TT 

an 

S 

Ss 





lelololololoelelelo) lejololololeloelelo) lelololelololelelto) 


Pleiotropy strength, ¢ 


0 02040608 10 02040608 10 02040608 1 


Fig 4. The evolution of pleiotropy promotes cooperation. Heatmaps show average trait values among the global 
population of cells (across all groups) at steady state in our model. Results are shown for 3 group sizes (increasing from 
top to bottom). Dotted line on each heatmap indicates separation between point where pleiotropy has no effect 
(control case ¢ = 0) to regions where it has a nonzero influence on the dynamics. Cooperation evolves for a wide range 
of parameters but is destabilised for longer group life spans, A, and larger group sizes, K, due to the emergence of 
cancerous cell lineages. This effect is strongest without pleiotropy (¢ = 0). When pleiotropy has an effect (¢>0), natural 
selection favours its evolution resulting in higher levels of cooperation, i.e., the blue intensity increases from left to 
right in the cooperation plot. As expected, pleiotropy is most favoured when it is more effective, i.e., the strength of 
pleiotropy, ¢, is higher. Note that the level of expression of the pleiotropic trait can be nonzero even when it has no 
effect due to stochastic effects. In these situations, however, the expression of the pleiotropic trait does not influence 
cooperative evolution. We show below that pleiotropy will also evolve with a cost, which greatly reduces this 
stochasticity (S11 Fig). Parameters: s° = s* = 0.95; K = 200; u = 0.0001; v = 0.01. The code required to generate this 
figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#. 
YjSBVurP2Uk. 


https://doi.org/10.1371/journal.pbio.3001626.g004 








Our analysis reveals a wide range of parameters where pleiotropy evolves and, in doing so, 
promotes cooperation over evolutionary time (Fig 4). In particular, when the effects of pleiot- 
ropy are absent in the model (¢ = 0), the levels of cooperation that evolve are often markedly 
decreased (Fig 4). These results includes cases where cooperation evolves to a lower level with- 
out pleiotropy but not to 0, which is also seen to a modest extent in the supplementary figures 
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of dos Santos and colleagues [29]. We also see many cases where cooperation is negligible 
without pleiotropy, but present with pleiotropy, predicting that pleiotropy can enable the evo- 
lution of cooperation in regions where kin selection alone does not maintain it (Fig 4, $8 Fig). 
Moreover, we do not see the evolution of pleiotropy in a control model of 2 private traits 
(Methods; S9 Fig). These 2 contrasts—removing pleiotropy and removing cooperation—dem- 
onstrate that cooperation can depend on the evolution of pleiotropy and vice versa. 

Our model predicts that pleiotropy is more important for cooperation as groups become 
larger and longer lived. This is because, if groups are small or short lived, there is less opportu- 
nity for cheater lineages to arise and interfere with group functioning (Fig 3). As expected, the 
strength of pleiotropy—how protective pleiotropy is against invasion of mutant cheater geno- 
types—is an important factor in determining when pleiotropy evolves (variation in the x- 
direction in Fig 4). Nevertheless, we find that pleiotropy evolves and increases cooperation 
even for low levels of protection (¢<0.5, blue plots, Fig 4). In some cases, pleiotropy evolves in 
the model when it has no phenotypic effect via genetic drift (Fig 4, ¢ = 0). As expected, this 
effect is strongest when group sizes are small and short lived because this is when natual selec- 
tion is also generally weakest. Below, we introduce a cost to pleiotropy that largely removes 
this effect and shows it does not contribute to our findings. 

Another assumption of potential importance is how detrimental the invasion of cheater lin- 
eages is to overall group function. If the invasion of even a few mutant cheaters can damage 
group function, then mechanisms that resist that invasion are likely to be more strongly favoured 
by between-group selection. To be conservative, we do not assume a high detrimental impact 
where only a few cheater cells are fatal for group functioning but instead consider a linear function 
where some multicellular function persists in the face of very large numbers of noncooperator 
cells. Specifically, recall that group function at age y is given by z;(y) = Z;,(y)zZ;,(y), a function that 
declines linearly with the invasion of noncooperative mutants (all else being equal). However, 
even without a high detrimental impact associated with a small numbers of cheaters, we find that 
pleiotropy is important for cooperation in groups that generate only 5,000 to 10,000 cells across 
their life span (i.e., 5 x 10° to 10* cell divisions in total). This prediction is borne out in Fig 4 
where the total number of cell divisions in the model corresponds roughly to group size multiplied 
by the number of generations in the life span (i.e., y-axis value multiplied by group size in Fig 4; 
see Methods). As discussed, these are only rough estimates as they depend on mutation rate and 
other assumptions that will vary between systems. However, given that many cellular groups 
undergo many more cell divisions than this estimate, these results predict that the evolution of 
pleiotropy has the potential for widespread impacts. 

Our model also allows us to follow the order in which the traits evolve in the population 
(Fig 5). Without pleiotropy (¢ = 0), cooperative genotypes can rapidly evolve, but so too do 
noncooperative genotypes (red lines in Fig 5A). As a result, these cheater lineages end up mak- 
ing up a substantial proportion of the population such that cooperation and group functioning 
are limited (Fig 5B, left panel). By contrast, when pleiotropy can influence the distribution of 
mutational effects, we see that it evolves extremely rapidly after the origin of cooperation itself. 
Indeed, for the higher strengths of pleiotropy, it evolves alongside the initial origin of coopera- 
tion in the population (Fig 5B, middle panel). 

Rather than being a secondary adaptation restricted to derived multicellular groups, there- 
fore, our model predicts that pleiotropy can be important at the origin of cooperation within 
multicellular groups. One can also assess the effects of pleiotropy in terms of trait-based muta- 
tional load: the average trait values of a group at age y relative to their trait value at birth (Fig 
5C). The mutational load is greatest for the cooperative trait, with groups tending to express 
reduced cooperation aged y than they do at their origin. However, this mutational load of 
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Fig 5. Long-term evolutionary dynamics of cooperation and pleiotropy. Dynamics are shown for the global population of cells over time, f, 
which encompasses many generations of cell groups. These dynamics encompass both within-group and between-group selection dynamics. We 
show 3 strengths of pleiotropy, ¢, to capture 3 scenarios where stronger pleiotropy is associated with more rapid and complete evolution of both 
pleiotropy and cooperation. (A) Changes in global genotype relative frequencies, x“(t). (B) Changes in the global average levels of cooperation, 
private trait expression, and pleiotropy, Z‘(t). (C) Average change in traits over a group lifetime, measured as the difference between the average 
trait values among groups aged y to those expected from their founding cell at birth. (D) Order in which genotypes invade in the global 
population, where a blue arrow indicates a gain of cooperation (via between-group selection), a grey arrow indicates a loss of cooperation (via 
within-group selection) and a black arrow indicates no change in cooperation. Parameters: s = s* = 0.95; K = 200; w = 0.0001; v = 0.01; 2 = 20. The 
code required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#. 
YjSBVurP2Uk. 


https://doi.org/10.1371/journal.pbio.3001626.g005 





cooperation is reduced with the evolution of stronger pleiotropy (Fig 5C). In this way, pleiot- 
ropy does not just act to reduce the emergence of cheaters in a given group (Fig 3), it acts across 
the whole population and can greatly increase the evolved level of cooperation (Figs 4 and 5). 


Pleiotropy evolves to stabilise cooperation across a wide range of 
conditions 


We have so far assumed that groups are formed from a single cell. While this is realistic for the 
majority of multicellular organisms, other cellular groups, particularly microbial groups, 
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commonly contain multiple genotypes that meet and mix. If large numbers of different geno- 
types meet and mix—and relatedness is close to O—the evolution of pleiotropy and indeed 
cooperation does not occur in our model. Under these conditions, there is no between-group 
genetic variation and the outcome of competition is determined solely by within group 
dynamics (S10 Fig). Here, so long as genotypes that lack the cooperative phenotype can arise 
at some point, they will take over and pleiotropy serves no function. However, relatedness can 
often be relatively high in microbial groups due to spatial structure, where a patchwork of 
groups form, each dominated by a single genotype [1]. We can study the effects of an interme- 
diate level of relatedness in our model by assuming groups are founded by 2 cells (chosen uni- 
formly at random from their parent group), such that there are now up to 32 different group 
genotypes in the population. This case has an important difference to the single-cell bottleneck 
case where cheater cells always start a new group alone with little chance of survival. With 2 
cells, cheater cells now have the chance of founding groups alongside cooperators that they 
can exploit, thus greatly improving their prospects. Despite the added complexity, we see again 
that the evolution of pleiotropy is often favoured and able to promote the evolution of coopera- 
tion as it evolves (S11 Fig). While the importance of pleiotropy in our model rests upon some 
relatedness between cells, therefore, it does not rest upon a single cell origin. 

Our conclusions are also robust to changing other assumptions and parameters. One key 
consideration is that there may be a cost to pleiotropy if, for example, the regulation of 1 trait 
is compromised by its linkage to another [41]. However, we find that pleiotropy still evolves if 
it carries such costs to a group’s functioning (S12 Fig), which is further testament to its ulti- 
mate importance for improving group function. Another important parameter is the benefit of 
cooperation (strength of group selection). Reducing the benefit of cooperation in our model 
reduces the scope for the evolution of cooperation but, importantly, where cooperation can 
evolve there are broad parameter ranges where pleiotropy evolves to increase cooperation (S13 
Fig). Notably, the evolution of pleiotropy is even seen when natural selection for cooperation 
is very weak, as may have occurred at the inception of multicellular life. 

Varying the relative probability of gain-of-function mutations has little impact on outcomes 
(S14 Fig). However, as expected, the baseline mutation rate is important. Increasing this muta- 
tion rate causes the more rapid breakdown of cooperation, which requires stronger pleiotropic 
effects for cooperation to be maintained. However, so long as strong pleiotropic links are possi- 
ble, we see that they rapidly evolve and again stabilise cooperation (S15 Fig). For reduced muta- 
tion rates, cheater lineages arise less often and so, even in the absence of pleiotropy, cooperation 
can be maintained more easily. All else being equal, therefore, pleiotropy will now only evolve in 
larger or longer-lived cellular groups. For example, halving the mutation rate (u = 0.00005) 
roughly doubles the number cell divisions where pleiotropy becomes critical for cooperation 
(compare Fig 4 with S16 Fig). However, even if we lower the mutation rate an order of magni- 
tude (u = 0.00001), we still observe the widespread evolution of pleiotropy in groups of only 
10,000 to 20,000 cell divisions (S16 and S17 Figs show this effect, with and without a cost to plei- 
otropy, respectively). In summary, we observe that the evolution of pleiotropy promotes cooper- 
ation for relatively small multicellular groups across a wide range of parameters. 


Pleiotropy is predicted to be important in a simple model of multicellular 
organisms 


Our modelling assumptions are most suited to multicellular groups of microbes and probably 
some of the ancestral organisms that gave rise to multicellularity in the algae, plants, animals, 
and fungi. With the evolution of derived multicellularity came many complexities, which our 
models do not capture. Central among these is the importance of the germ and soma 
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Fig 6. The evolution of pleiotropy in multicellular groups with a germ line. (A) Cheater lineages that threaten 
cooperation evolve in multicellular species with a germ line, as well as in simpler multicellular groups. In multicellular 
organisms, these lineages manifest themselves as tumours and cancers that emerge by mutation. From left to right, 
bacteria, simple and complex animals. Bacteria: Spontaneous /JasR mutant emerge during evolution of the pathogen 
Pseudomonas aeruginosa. The mutants do not contribute to the production of public goods molecules that their wild- 
type counterparts do, which enables them to gain a short-term competitive advantage (Image credit: Sheyda Azimi and 
Steve Diggle). Simple animal: Spontaneous tumour formation in the basal metazoan Hydra, a tumour-bearing Hydra 
oligactis polyp (right) is shown next to a healthy animal (left). Credit: Alexander Klimovich, Kiel University. Complex 
animal: Cancer is widespread in long-lived animals, especially zoo animals where other selection pressures are 
minimised. Nanuq, a 29-year-old polar bear, died from liver cancer in The Columbus Zoo, Ohio, in 2017 (image shows 
a different polar bear; credit: Alan D. Wilson). (B) To study the impact of pleiotropy on multicellular groups with a 
germ line, we varied y, a measure of the strength of the transmission of the germline, modelled as the likelihood that a 
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group propagates the genotype of its founding cell at the age at which it reproduces versus a cell selected at random. 
Heatmaps show average trait values among the global population of cells (across all groups) at steady state in our 
model. Results are shown for 3 germ line strength parameters (increasing the strength of the germline transmission 
between ancestor and descendent groups from top to bottom). Increasing the strength of the germline has a positive 
effect on the evolution of cooperation, but pleiotropy continues to be favoured in regions of the parameter space in 
which cooperation is vulnerable to breakdown. The dotted line marks the boundary between pleiotropy having no 
effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: s“ = s* = 0.95; K = 200; u 
= 0.0001; v = 0.01; ¢ = 0.02. The code required to generate this figure can be found at https://github.com/euler-mab/ 
pleiotropy and https://zenodo.org/record/6367788#.YjSBV urP2Uk. 


https://doi.org/10.1371/journal.pbio.3001626.g006 





separation for modern multicellular life. We, therefore, next ask whether this biology influ- 
ences the evolution of pleiotropy and multicellular cooperation. 

Here, we assume that multicellular groups can influence which genotype they transmit dur- 
ing reproduction via a germ line that undergoes fewer cell divisions, and mutations, than the 
soma. Specifically, with probability y, an ancestral group seeds a new group with a cell of the 
same genotype to its founding cell, and with probability 1—y, a cell is chosen uniformly at ran- 
dom from the group to seed a new descendant (as before). As discussed above, we also assume 
a cost to pleiotropy to be conservative. The model predicts, as expected, that a germ line is gen- 
erally benefitial for the evolution of cooperation, as it reduces the chance that noncooperative 
genotypes will start new groups. Nevertheless, with a germ line, we still see the widespread evo- 
lution of pleiotropy (Fig 6). 

A second important characteristic of derived multicellular organisms, like humans, is the 
need to reach a certain age before reproduction is possible. Up until now, reproduction among 
our multicellular groups has been age independent. We therefore introduce another parame- 
ter, a, which determines the fraction of the expected life span that groups must reach before 
they can reproduce. For example, if the expected life span is 2 = 50 and a = 0.5, then groups 
can only reproduce after age y = 25. On its own, adding this requirement reduces the levels of 
cooperation in the population, because it means that organisms reproduce when they have the 
highest levels of cheater mutations (S18 Fig). However, the negative effects of late life repro- 
duction on cooperation are reduced if we assume that these organisms also have a germ line 
(Fig 7). Moreover, this effect rests upon the ability of an organism to evolve effective pleiotro- 
pic constraints (¢ > 0). We find, therefore, that organisms with a germ line and delayed repro- 
duction will evolve pleiotropy to protect against cheater lineages and ensure high-level 
functioning when they reproduce (Fig 7C). 


Discussion 


Our models predict that pleiotropy can greatly increase the levels of cooperation within multi- 
cellular groups (Fig 4 and S11-S18 Figs). Natural selection favours pleiotropy because it 
reduces the probability that cheater lineages arise within a cellular group or multicellular 
organism, which would otherwise damage functioning (Fig 3). We expect this process to be 
important whenever groups become large or old enough for the emergence and spread of 
cheater lineages to cause significant harm to multicellular function (Fig 4). These conditions 
are expected in relatively small and simple multicellular groups. We find pleiotropy can evolve 
to suppress cheater lineages in groups that undergo as few as 10° cell divisions. 

This 10* cell divisions estimate depends on parameters, most notably mutation rate. Muta- 
tion rates vary widely, but to give an idea of natural rates, a simple animal like Hydra vulgaris 
would be expected to have between roughly 1 and 100 mutations per genome per cell division, 
based upon its approximately 1 Gbp genome and recent estimates of somatic mutation rates 
taken from mammals (10°-10” mutations per basepair per cell division) [42,43]. If we divide 
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Fig 7. The evolution of pleiotropy in multicellular groups with a germ line and delayed reproductive maturity. 
We varied a, the fraction of the expected life span A groups must have aged to before they can reproduce. We also 
assumed a germ line, y = 1, and that the evolution of pleiotropy carries a cost. Heatmaps show average trait values 
among the global population of cells (across all groups) at steady state in our model. Results are shown for 3 
reproductive maturity parameters (increasing the age at which maturity is reached from top to bottom). Increasing the 
age of reproductive maturity favours the evolution of pleiotropy. The dotted line marks the boundary between 
pleiotropy having no effect (control case) and pleiotropy having an effect on the outcome of mutations. Parameters: s“ 
= s* = 0.95; K = 200; w = 0.0001; v = 0.01; € = 0.02; y = 1. The code required to generate this figure can be found at 
https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2UKk. 


https://doi.org/10.1371/journal.pbio.3001626.g007 





these mutation rates by the number of genes in H. vulgaris (20,000), this gives rates between 
one-half and 50 times our standard value of y (0.00005 to 0.005). For a trait affected by a single 
gene, this range may still overestimate the rate of loss-of-function mutations, because not all 
mutations in a gene will hit a coding region nor lead to a loss of function if they do. However, 
the traits we are interested in here—for example, regulated cell proliferation—are often 
affected by large networks of genes, such that mutations in many genes have the potential to 
influence the trait, for example, drive unregulated cell proliferation. Our standard value, there- 
fore, appears to be reasonable, but in species with lower mutations rates, a larger number of 
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cell divisions are required before pleiotropy is needed to maintain cooperation (S16-S17 Figs). 
However, even if 10* cell divisions is an underestimate by 1 or even several orders of magni- 
tude, many multicellular groups will remain above the threshold at which pleiotropy is pre- 
dicted to be important. There are an estimated 10° cells in hydra [44], 10° neurons in the 
Drosophila brain [45], 10°° cells in a bacterial colony [46], and 10° cells in humans [47], where 
these large standing cell numbers can still greatly underestimate total cell divisions because of 
cell turnover [30]. If relatively few cell divisions are needed for the evolution of pleiotropy, this 
suggests that pleiotropy may have had an important role in cheater suppression at the origins 
of multicellularity. Further consistent with this hypothesis, the evolutionary dynamics in the 
model predict that pleiotropy will evolve rapidly, close to simulatenously, with the first evolu- 
tion of cooperation (Fig 5). 

We have identied a scenario where pleiotropy evolves because of its positive effects on 
cooperation, but pleiotropy can also arise for other reasons [25,35,39,48]. Pleiotropy is wide- 
spread for all types of traits—cooperative and noncooperative—and can readily arise, for 
example, as a way to coregulate multiple traits [39,48]. Such regulatory structures might help 
to stabilise a cooperative trait or help cooperation to first arise when there is positive fre- 
quency-dependent selection [49]. A challenge for the future is to distinguish between examples 
where the effects of pleiotropy on cooperation arose independently of cooperation, and cases 
where cooperation drove the evolution of pleiotropy. Our within-group model underlines 
that, no matter why pleiotropy first evolves, it has the potential to promote cooperation (Fig 
3). Moreover, when individuals are related, this can generate between-group selection for plei- 
otropy and the stabilisation of cooperation in the long term (Figs 4 and 5). 

There are other evolutionary processes than those we have modelled here that may lead to 
cooperation being associated with pleiotropy. Whenever cooperation is linked to the ability of 
a species to compete and persist in an ecosystem [50], for example, species-level selection may 
enrich for species that have evolved pleiotropic links that promote cooperation over those that 
do not [51]. Another potential route to pleiotropy is when 1 individual evolves to enforce 
cooperation in another. For enforcement to be effective, the enforcer needs to find a way to 
constrain the recipient in some way from escaping the enforcement. This need can result, 
therefore, in a pleiotropic link in a recipient where the cessation of cooperation is linked to a 
personal cost from enforcement [25,28]. This case is notable in that it does not rest on there 
also being positive relatedness between interacting individuals (cf., S10 Fig). 

Enforcement may be important in microbial examples of pleiotropy, including the case of 
Dictyostelium discoideum [24] where enforcement of 1 cell type by another appears to be 
important for cooperation [52]. However, our modelling also suggests that pleiotropy may 
evolve in microbes to limit the spread of cheater lineages. Mutation rates can be relatively low 
in microbes [53], and the arrival and mixing of multiple genotypes within a microbial group 
has the potential to both undermine the evolution of cooperation and the potential for pleiot- 
ropy with it [29]. However, cooperative traits are known to be widespread in microbes, partic- 
ularly in the well-studied bacteria [54]. Moreover, groups of bacteria reach vast numbers and 
the evolution of noncooperator lineages has been observed [21,55,56]. Whenever these line- 
ages cause significant harm to cooperative function, therefore, there is the potential for pleiot- 
ropy to evolve to prevent their emergence. 

We predict that the evolution of pleiotropy will be particularly important for multicellular 
organisms. Development from a single cell (Fig 4 versus S11 Fig), combined with relatively 
high somatic mutation rates (S15-S17 Figs) and the potential for a very large numbers of cell 
divisions (Fig 4 and S16 and S17 Figs), are all predicted to favour pleiotropy as a means to pro- 
mote cooperation. Consistent with this prediction, many multicellular organisms suffer from 
cancer, which has long been viewed as analogous to the evolution of cheater lineages within 
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cellular groups [29,31,57-61]. Moreover, the study of cancer has identified powerful mecha- 
nisms to supress potential cheater lineages, which function by linking oncogenic traits (loss of 
cooperation) to secondary traits that limit viability (loss of private trait) [31,38,62,63]. For 
example, loss of function mutations in the retinoblastoma tumour suppressor gene are linked 
to the activation of programmed cell death (apoptosis) [36,37], while activation of ras genes— 
key regulators of cell proliferation—are associated with cellular senescence and cell cycle arrest 
[64]. Other mechanisms push cells with oncogenic mutations to differentiate from a prolifer- 
ative cell type (stem cell) to one that does not divide [63]. The links between cooperation and 
cell viability can be cell intrinsic but they also occur via the actions of immune cells and othe 
cell types [62,65], which can promote apoptosis in a mutant cell [66,67] or even physically 
force it out of a proliferative tissue [68]. The biology of multicellularity, therefore, appears to 
be characterised by a large number of pleiotropic connections between abberant cell prolifera- 
tion that threatens cooperation and the loss of traits that cells need to survive and divide. 

The effects of pleiotropy in our models is testament to the importance of considering 
molecular mechanisms in the study of cooperation, particularly in cellular systems [31]. In 
microbes, for example, mechanisms such as prudent regulation of cooperative traits, quorum 
sensing, and green beard genes can all serve to promote cooperation [69-71]. High relatedness 
is also clearly important for cooperation across a range of cooperative systems, and is particu- 
larly important for cooperation in cellular groups. Experiments have demonstrated how relat- 
edness promotes cooperation in microbes, and development from a single cell was likely to 
have been central to the evolution of complex multicellularity [14,15]. However, the large 
numbers of cell divisions that occur in many cellular groups means that single-cell ancestry is 
far from sufficient to maintain cooperation. Indeed, our models predict that even relatively 
small and simple cellular groups will suffer if cheater lineages are left unchecked. This suggests 
that, in addition to the single-cell bottleneck, the evolution of pleiotropic constraints may have 
been important for the origins of multicellularity and the subsequent major transition in evo- 
lution that led to the complex multicellular life that we see today. 


Methods 
A general model of multilevel selection using an age-structured model 


Capturing the full dynamics of a multilevel selection process is difficult because the potential 
for selection within and between groups creates an enormous space of possible outcomes. For 
these reasons, many early models of multilevel selection, such as the haystack model [40], 
attempted to simplify the problem by neglecting to model group-level events explicitly. While 
this simplifies the model, the approach misses the fact that group-level events and individual- 
level events can happen on different timescales, something that is particularly important for 
capturing the biology of most multicellular groups. 

To capture the 2 levels of selection explicitly, we study the dynamics of natural selection in 
an age-structured population of multicellular groups and in subpopulations of cells within 
those groups. Changes in the distribution of groups of different ages and types evolves with 
time ft and age y. Changes in the distribution of cells of different types within each group 
evolves with age y. Using a similar approach to [72], we use a system of PDEs to model changes 
in the age-structured population of different types of groups over time. We model changes in 
abundances of different cell types within the different group types as they age using stochastic 
simulations. These simulations then allow us to characterise and describe how groups of differ- 
ent types differ in their development as they age. Within-group differences in development 
further provide the basis for differences in reproductive success between groups founded by 
different types of cells. All the code used to perform our numerical analyses is open source and 
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can be found in GitHub (https://github.com/euler-mab/pleiotropy) and Zenodo (https:// 
zenodo.org/record/6367788#.YjSBVurP2Uk). 

We begin by deriving and describing the structure of our model in general terms, before 
turning to the specifics of how we use it to study the evolution of cooperation and pleiotropy. 


Cell and group types 


We assume there is finite set of cell types, i,j€I, and a finite set of group types, k,IcK. We refer 
to cells of type i as i-type cells and multicellular groups of type k as k-type groups. To distin- 
guish between other group-level and cell-level variables in our model, we use superscripts g 
and c, respectively. The type of a group is assumed to be defined by the type(s) of its founding 
cells. In our standard model, we assume that a group is founded by a single cell. In this case, 
the set of cell genotypes and group genotypes is the same (i.e., I= K). In a later model, we 
assume that a group is founded by 2 cells. In this case, the set of group types contains all unor- 
dered pairs of cell types, and the index notation becomes more cumbersome, although all the 
same principles for the model hold true. We therefore focus our description here on the sim- 
pler model. Let us now describe how the abundances of cells and groups changes over time. 


Between-group dynamics 


We begin with the group population. In the limit as the population of groups gets large, we 
assume that relative density of k-type groups aged y can be modelled as a continuous quantity, 
n(t,y) € [0, 1]. In the absence of births and deaths of whole groups, all groups in the popula- 
tion simply age. What this means is that if there were a population of 10 groups at time t aged 
y = 1, then at time t+10, there will be 0 groups aged y = 1, and 10 groups aged y = 11. 

Consequently, the relative density of k-type groups aged y at time ¢ changes over time 
according to a system of PDEs satisfying a conservation law of the form 

Ony, : 


Tet) + GE ley) =O (1 


where there is no change in the overall density of groups in the population. This law can be 
derived as follows. First, note that within a particular age range [y,, yz], the total abundance of 
k-type groups aged y at time t is given by 


[ ileney (2) 


If we assume that there are no births or deaths of groups within the age range, then the 
abundance of individuals in the age range [y), y2] can only change because of a process of age- 
ing. Groups of a younger age may enter this age range at the lower age y,, and those within the 
age range may get older than the upper age bound y,. If we suppose that groups age at a con- 
stant rate v, then the rate of change of k-type groups aged y at time f is just vjj(y, t). The rate 
of change of the total abundance of individuals at time ¢ in the age range [y,, y2] is given by 


v2 


d 
dt ny(t, y)dy = v(ni(t, 1) _ ni(t,¥)) (3) 
pal 


where it is equal to the flux of the ageing population over the boundaries of the age range. Let 
us now integrate both sides of this equation to get an expression for the abundance of k-type 
groups in [y;, yz] at time t)>t, in terms of the abundance of k-type groups aged y at time f, 
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and the total flux at each boundary during this time period. We have 


i ” (tty) — nity y) dy = —v / “(nf(t,y,) — nf(t,y,))de. (4) 


M1 


Assuming n;(t, y) is differentiable, then we can use the fundamental theorem of calculus to 
rewrite the integrands on both sides as 


df? 29 
riley) — Mle =5 fmt yar fs me,ydde, (5) 
and 
HA (t,y») — n(ty,) = 4 ‘i ” a(t, y)dy = / "2 aklt,y)dy (6) 
kha S2 kM Po Y1 By a Oe o) 


respectively, where we have used Leibniz integral rule to take the derivative operators inside 
the integrals on the right-hand side. Substituting back in to Eq (4) and rearranging, we have 


a S Rey iy? hE Vaya (7) 
is : Ot k y oy k Wy y= : 


Since this integral is 0 for any arbitrary age range and any time interval, we must conclude 
that the integrand itself is exactly 0: 


0 & 0 & aa, 
py teltsy) + Va ee = 0. (8) 


We can always normalise the rate of ageing to a constant v = 1, meaning that this equation 
simplifies to the conservation law given by Eq (1) we started with above. 

These equations are similar to a PDE derived by Burt Simon (see [72]), which was also used 
to study group selection. An important distinction is that our approach allows one to study 
group-selection using a system of n PDEs, rather than 1 PDE evolving on an n-dimensional 
surface. The system of PDEs can simplify the numerical problem of solving group-selection 
dynamics dramatically but is generally only tractable for biological systems in which groups 
can be assigned to categories in some simple way. Here, for example, we categorise groups by 
the genotype of their founding cell. If groups were formed from multiple founding cells, then 
the number of group types would increase dramatically, making the numerical solution of 
these equations much more difficult. 

To incorporate assumptions about the birth and death of k-type groups in our PDE model, 
along with their production by other /-type groups due to mutation, we can simply extend Eq 
(1) to include source and sink terms. The production of k-type groups aged y due to the repro- 
duction of all other /-type groups can be represented by a term 


ee) 


So ni(t.y OF (ty hy) ay’ 


0 le{I} 


where bi(t, y’) € [0, co} is the rate at which /-type groups aged y’ reproduce at time ¢, and 
hi.(y’,y) € [0,1] is the conditional transition probability that a k-type group aged y is pro- 
duced, given that an /-type group aged y’ just reproduced (in practice, we assume the age of a 
newly produced group is always y = 0). The death of groups can simply be modelled by a term 


—ny(t, y)di(t,y) 
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where dj(t, y) € [0, co} is simply the rate at which k-type group aged y die at time ft. Incorpo- 
rating these terms into Eq (1) gives 
& 


vs Om 


& foe) 
ay) -| Soren by OEY Mi yay — mit y)ai(t.y)->*(ty). (9) 
Ot 0 e{T} Oy 


The extent to which there is between-group competition in our model therefore depends 
on variation between different types of groups in birth rates, bi(t, y), death rates, di(t, y), and 
mutation transition probabilities, /j,(y, y/). These functions all depend on the within-group 
dynamics within each type of group. 


Within-group dynamics 

Now let us describe what happens within groups as they age. We assume all groups are 
founded by a small number of cells (1 cell in most of our analyses) and undergo logistic growth 
as they age. The abundance of i-type cells within a k-type group aged y is denoted nj,(y). We 
assume that external forces have no impact on the change in abundance of cells of different 
types within a group. Thus, the rate of change in the abundance of i-type cells within a k-type 
group aged y depends only on age y, and not on external time t. 

Changes in the abundance of different cell types within a particular k-type group are sto- 
chastic, but we assume that all groups of type k generally develop in the same way as they age, 
irrespective of the environment they were born into. In practice, we therefore characterise the 
expected development of a k-type group by calculating the average behaviour of 10,000 repli- 
cates of the stochastic dynamics of cells within the group. We use stochastic simulations rather 
than deterministic solutions to characterise the within-group dynamics because we are inter- 
ested in the different rates at which mutant cell lineages invade within a finite subpopulation. 
If we used ordinary differential equations to characterise the same behaviour, cheater mutants 
would emerge and spread within each group deterministically very early on in the lifetime of 
groups due to the assumption of infinite population sizes. In contrast, our averaging approach 
shows that in finite systems, mutant cell lineages invade at vastly different rates in different 
types of groups. This is crucial for showing why pleiotropy is important in group selection. 
Groups founded by pleiotropic cooperators can outcompete groups founded by nonpleiotro- 
pic cooperators because pleiotropy slows the rate of invasion of mutant cell cheaters, making 
groups more competitive. 

The stochastic dynamics within groups are characterised by a birth-death process with 
mutation. The birth rate of i-type cells in a k-type group aged y is given by bj, € [0, 00}. The 
death rate of i-type cells in a k-type group aged y is given by dj, € [0, 00}. Finally, the condi- 
tional transition probability that an i-type cell is produced, given that a j-type cell just repro- 
duced (akin to a mutation transition probability), is denoted hi, € [0, 1]. With these rates 
defined, we can simulate the stochastic evolutionary dynamics of any population. We generate 
sample paths using Gillespie’s Direct method [73]. 


Solving the equations 


So far we have described how selection and mutation acts on groups and cells within groups. 
How are these equations linked? The answer is that the group birth rate, bi(t, y’), death rate, 
di(t, y’), and transition probability, hi, (y’, y), all depend on the outcome of the within-group 
cell dynamics. Thus, to solve the system of PDEs describing the dynamics of selection and 
mutation between groups, we first need to resolve the dynamics of selection and mutation 
between cells within each type of group. 
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We resolve the within-group dynamics of cells using custom-built stochastic simulations 
written in Matlab (R2019b). We then numerically solve the between-group dynamics given by 
Eq (9) using a custom-built finite-volume PDE solver written in Matlab (R2019b). We use a 
“superbee” flux limiter to control the rate of change of group relative density in regions where 
the group density function is not smooth [74]. 


Modelling the evolution of pleiotropy and cooperation 


Now that we have described a general model of multilevel selection in an age-structured popu- 
lation, we turn to describing how we use that model to study the evolution of cooperation and 
pleiotropy. Specifically, we focus on describing the birth, death, and mutation rates for cells 
and multicellular groups in our model. For reference, all of the parameters used in the model 
and their default values and ranges are shown in Table 1. 


Cells 


Cells in our model express 3 phenotypic traits of interest, which can either be in an active or 
inactive state. There are therefore 8 possible cell genotypes in our model (Fig 2B). We consider 
a cooperative trait z° € {0,1}, which is beneficial for the reproductive success of multicellular 
groups as a whole, but costly for the reproductive success of cells within the lifetime of a group, 
a private trait, z¢ € {0, 1}, which is beneficial for the reproductive success of a cell expressing 
it within a group and beneficial for multicellular groups as a whole, and a “pleiotropy” trait, 
Ze {0, 1}, which influences the types of cell mutations that can occur. Ordering cell traits as 


row vectors, [z/,, z,, Z;], the genotype-phenotype map for cells can be represented by a matrix 






































0 0 0 
0 0 1 
0 1 0 
01 1 
z= . (10) 
1 0 0 
1 0 1 
1 1 0 
1 11 
Table 1. List of parameters used in the model. 
Parameter [= default Description 
value] 
s€ = 0.95 Strength of within-group selection. 
s¥ = 0.95 Strength of between-group selection. 
A= 20 Expected life span of groups (units of time). Note: The expected life span of a cell is 
roughly 1 unit of time. 
K= 200 Carrying capacity of groups (size at maturity). 
u = 0.0001 Loss-of-function mutation rate of traits. 
v=0.01 Relative rate of gain-of-function mutations to loss-of-function mutations. 
¢ =0.0 Strength of pleiotropy. 
€=0.0 Cost of pleiotropy (% reduction of group function). 
a=0.0 Age of reproductive maturity. 
y = 0.0 Strength of germ line transmission. 
https://doi.org/10.1371/journal.pbio.3001626.t001 
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Each of the 8 rows corresponds to a genotype, and each of the 3 columns corresponds to a 
trait, ordered public (cooperative), private, and pleiotropy. Cells with different phenotypes 
may vary in their ability to survive and reproduce within groups. The within-group birth and 
death rates of i-type cells in a k-type group aged y are given by 


Cyl 
1. — sz, 


ae eel a 
i) 1 _ sZ,(y) ’ 


(11) 


and 


(1— sei) (Ney) — 1) 
(=F) K 





di(y) = (12) 


respectively, where s‘ is the strength of selection on cell traits, z°,(y) and Z‘,(y) are the average 
expression of the cooperative and private traits in a k-type group aged y, respectively, and 
N;,(y) is the size of a k-type group aged y. The component N;(y) — 1 in the death-rate prevents 
a group from dying because of stochastic extinction of all its cells, because its death rate is 
d;,(y) = 0 when the size of the group is Nf(y) = 1. 

Note that the expression of the cooperative trait z‘, places a cost on i-type cell division rela- 
tive to other cells in the group, the expression of the private trait, z‘,, gives i-type cells a survival 
advantage relative to other cells in the group, and z/; has no impact on cell birth rates or death 
rates at all. We are assuming here that the cooperative trait positively impacts the birth rate of 
cells and the loss of the private trait positively impacts the death rate of cells. These assump- 
tions are based on the observations that extracellular growth factors and enzymes in cellular 
groups (cooperative traits) often promote growth [1], while cells lacking a functional metabo- 
lism (a private trait) often die [75]. 

We generally expect within-group selection to disfavour expression of the cooperative trait 
z‘, favour the expression of the private noncooperative trait z‘,, and be neutral with respect to 
the pleiotropic trait z,, Similar to dos Santos and colleagues, we assume that selection is of 
equal strength with respect to the cooperative and private traits because we did not want to 
introduce biases by privileging 1 trait over another, but rather focus on the importance of plei- 
otropy [29]. In a variant of our model, we switch the cooperative trait to another private trait, 
z, € {0,1}. To do this, we replace our usual birth rate function with 


ic C py 
Ll-s+sZ, 


b.. S —_ ——— 
ui) 1—s°+s%,(y)’ 


(13) 


where the expression of the private trait z°, on i-type cell division is positive relative to nonex- 
pressing cells in the group. 

Cell mutations can occur during cell division. The conditional transition probability that an 
i-type ancestral cell produces a j-type descendant cell, given that it has reproduced is given by 


hi, = opi, + (1 — 6) 45 (14) 


where $¢[0, 1] is a parameter scaling the strength of pleiotropy, and p‘, and qi, are elements of 
2 different mutation matrices, the first of which represents a model of pleiotropic mutations, 
and the second of which represents a null model where pleiotropy has no effect. The parameter 
¢ thus scales the likelihood that mutation rates are sampled from a mutation matrix in which 
pleioptropy has an effect versus a mutation matrix in which pleiotropy is absent. When ¢ is rel- 
atively small, mutations to traits are almost always independent events, and pleiotropy is 
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therefore relatively weak or absent, but when ¢ is close to unity, mutations are often noninde- 
pendent events, where mutation of 1 trait influences the state of another trait. 

The specifics of our asssumptions about the effects of pleiotropy are as follows. Under our 
null model, we have a mutation matrix Q° = [qj], where we simply assume that mutations 


affect each trait independently. We assume that loss-of-function mutations occur with rate 
and that gain-of-function mutations occur with rate vu, where v<1. The mutation matrix for 
this model is then given by 



































(1 — vu)" (1 — vu) vq (1 — vu)"vn (1 — vu)? (1 — va)’va (1— vu) ee (1— vu)" we vie 
(1—va)'n (L—vu)(—w) (Ava) (A va)(1—p)ye (vu) (1 v)(1— wu vie (1- we 
(l1—vu)’e = (Lava)? (1 va)'(1—n) (L—va)(1— eye (1 vg? ve (L—p)(1—va)ye (1 wv? 
: — vy) pe vu Ku vu HE Vil mu) ve — pve — pve — pn)’ ve 
O= (l—vu)w A -vw)d-wWe d-w)dl-we (1-4) (1-4) (1 — pve (1-4) (15) 
(1 vu) (1 vp)vee (1 — vu) ve? vie (1—va)(1—w) (1—vw)(1- ave (1 ve)(1— ayn (1 we 
—vp)e (1 ve). — pe vee — Wve — ve)(1— we = vp)(1 =p)? — pve — yu) vu 
(1—vu)w’ (1 — va)(1 — 4) (1 = nye (L—-vu)(1—- we (1—va)(L— 4) (1 = wv (1-4) 
— vy) pe ve —ve)(1— p)e — nye —vu)(1 = w)u — wy? —vu)(1—p)’  (1—)’ve 
(1— vu) (1 — vu)(1 — p) (1 — pve (1— vu)(1 — H) (1-4) (1— vu)(1 — 4) 
Lo (- we (l- we (1—4)'u (l- we (l—)'h (1—)*n (1-4) J 
Under our model of pleiotropy, we have a mutation matrix P* = [pj], where we assume that 
a loss-of-function mutation to either z‘ or z° drives a corresponding change in the other trait 
(i.e., that loss-of-function mutations have pleiotropic effects). However, we assume that pleiot- 
ropy is only active in the cell genotype in which all 3 traits are active (i.e., pleiotropic coopera- 
tors with phenotype z,, z;,z, = 1). The mutation matrix for this model is then given by 
[ d—va) (vm) ve (1 —vu)’ve (1 = vpy'e (1 —vp)’vm (1 = vay’ (1 = ype vie 
(lve) (l— vey) (vp) (Lv) pve (Lave)? (1 = va) — wv vw (l= wn 
(L—vu)'e  (L— vey (1 vp)(L= pn) (A ve)(L— pve (1 vy vue (1— w)(. — va) (= pve 
e (L-va)ye (-ww)(l-we (-w)d-we (= ve)(l— 2) vee (1- pyr (l- py? (1 = p) vu (a6) 
(1—vy)'a (1 —vp)vp? (1 = vp)ve? vue (1—va)'(L— ph) (1—ve)(1—p)ve (1 vp)(1 — pve (= pve 
(l= vy)e? (1 ve) pe vee (1-p)vu (l-vu)(1— pe (L—ve)(1 = pn)? (1— py? (1 = p)*vu 
(1 — vu)? ve (1 = vu)(1 = 2) (l- py? (1— vu). — wu (1— pyr (1 = vu)(1 =p)? (1 = p)*vu 
2p ye 0 0 0 0 0 (l-n)'n 1—(2u— yp?) — (1— pw)" 





where the only difference between P* and Q‘ is in the final row. 


Multicellular groups 


Groups in our model express their own phenotypic traits of interest that emerge from interac- 
tions of cells within a group. We usually assume that groups are founded by a single cell, in 
which case there are 8 group genotypes in our model. In a variant of our model, however, we 
assume that groups are founded by 2 cells. In this scenario, there are 8(8—1)/2+8 = 36 group 
genotypes. We assume that both the cooperative and private traits of cells are essential compo- 
nents of the ability of the groups to function as a cohesive whole in order to maximise survival 
and reproduction of the group. The functionality of k-type groups aged y is denoted z;,(y) € 


(0, 1] and defined as 
(17) 


ZY) = Zn (V)Zru(Y)> 


where there is a multiplicative effect of the average expression of the cooperative trait, Z°,(y), 
and average expression of the private trait, Z°,(y), among cells within the group at age y. It is 
important to note that whereas cell-level traits are all fixed within their lifetime, the 
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functionality trait z;,(y) is expected to change throughout the lifetime of a group as it ages. We 
usually anticipate that within-group mutation-selection dynamics among cells will degrade 
group functionality as a function of age. In 1 variant of our model, we add in a cost of pleiot- 
ropy on the functioning of groups by defining 


ZY) = Zu (V) Ze) (1 — CZu(¥)), (18) 


where the component 1 — cam (y) allows us to explore what happens to the evolution of pleiot- 
ropy under the assumption that it carries a cost, ¢. 

Whereas the expected growth rate of cells is implicitly defined as 1 cell division per unit 
time, we assume that groups have an expected life span 2>1, which gives an expected group 
reproduction rate of p = 1/A. The larger A is, the more cell divisions a group will experience 
within its life span. We also consider that group reproduction might depend on reproductive 
maturity, achieved at a fraction of the expected life span, a€[0, 1]. For example, if a = 0.5 and 
A = 50, then the age of group reproductive maturity will be 25 and groups will only be able to 
reproduce after that point. The reproductive maturity status of a k-type group aged y is 
denoted z,(y) € {0, 1} and defined by 


oe eae (19) 


The functionality trait and reproductive maturity trait influence the birth and death rates of 
groups within the group population. Specifically, the between-group birth and death rates of 
k-type groups aged y at time t are given by 


p(l ~ st + s*zn(y)) 





BUY) =" ary weee(H)) (20) 
and 
: _ p(l — s8za(y)) 
di(t,y) = Tas)” (21) 


respectively, where s* is the strength of selection on group traits, Z*(t) is the fraction of multi- 
cellular groups in the group population that have reached reproductive maturity at time t, and 
Z;(t) is the average group functionality in the population at time f. Note that groups with a 
higher level of functionality gain a survival advantage over groups with a lower level of func- 
tionality. Thus, between-group selection favours maximal expression of both cooperative and 
private traits among cells. 

An ancestral group will sometimes produce a descendant group with a different set of 
founding cell genotypes than it started life with. By default we assume that a descendant group 
is formed by choosing a founding cell at random from the ancestral group at the time at which 
it reproduced. Under these circumstances, the conditional transition probability that a k-type 
ancestral group produces an I-type descendant group, given that it has reproduced, is given by 


haly, 0) = Sei) Min (22) 


where x;,(y) = nj,(v)/N{(y) is the relative frequency of I-type cells in a k-type group aged y 
and hy, is the conditional transition probability that a k-type cell mutates to an /-type cell dur- 
ing the production of the founding cell of a new group. Of course, Hj,(y, 0) is nonzero and so 
the most likely scenario is that a cell does not mutate. Note that nj,(y) and N;(y) are the 
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absolute abundance of I-type cells in a k-type group and the size of a k-type group aged y, 
respectively. 

In our default scenario, the production of mutant descendant groups tends to become more 
likely as a group ages, because mutation-selection dynamics among cells degrade its clonality. 
In a later extension of the model, we consider the alternative possibility that groups might 
have a germ line. To model this scenario, we assume that with a probability y, groups are 
founded by a cell with the same genotype as the ancestral group (except for rare mutations in 
germ line cells), and with probability 1—y, they are founded by cells chosen at random with 
probability proportional to their frequency, as in our default scenario. The conditional transi- 
tion probability that a k-type ancestral group produces an /-type descendant group, given that 
it has reproduced, in this scenario is given by 


hi(y, 0) = phy + (1-1) > pes”) Fi (23) 
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S1 Text. Revisiting the analysis of dos Santos and colleagues [29]. This Supporting informa- 
tion goes through a previous study on the evolution of pleiotropy and concludes that a differ- 
ent modelling approach is needed. 
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$1 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 

i= 7. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype g,, which actively expresses cooperation, z‘ = 1, a private trait, z¢ = 1, but no pleiot- 
ropy; Z, = 0. Growth of the group as its age, y, increases, is logistic, with a carrying capacity 

K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ¢. The 
vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% 
of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. 
(A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative frequencies, 
x(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, 
Z‘(y). (D) Differential fitness effects of loss-of-function mutations within the group. For exam- 
ple, the blue bar indicates the fitness effect of a loss-of-function mutation in the cooperative 
trait. Parameters: s° = 0.95; K = 200; u = 0.0001; v = 0.01. The code required to generate this fig- 
ure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/ 
6367788#.YjSBVurP2Uk. 

(DOCX) 





$2 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 
i= 6. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype ge, which actively expresses cooperation, zi, = 1, pleiotropy, z; = 1, but no private trait, 
z° = 0. Growth of the group as its age, y, increases, is logistic, with a carrying capacity K = 200. 
Dynamics are shown from left to right for 3 different strengths of pleiotropy, ¢. The vertical 
dashed line in (A-C) represents the point at which mutant cell lineages make up 25% of the 
group. Note that the strength of pleiotropy has no effect on the within-group dynamics. (A) 
Changes in genotype abundances, n‘(y). (B) Changes in genotype relative frequencies, x“(y). (C) 
Changes in the average levels of cooperation, private trait expression, and pleiotropy, Z‘(y). (D) 
Differential fitness effects of loss-of-function mutations within the group. Parameters: s° = 0.95; 
K = 200; u = 0.0001; v = 0.01. The code required to generate this figure can be found at https:// 
github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 
(DOCX) 
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$3 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 

i= 5. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype gs, which actively expresses cooperation, z° = 1, but no private trait, z° = 0, and no 
pleiotropy, z; = 0. Growth of the group as its age, y, increases, is logistic, with a carrying 
capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, 
¢. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 
25% of the group. Note that the strength of pleiotropy has no effect on the within-group 
dynamics. (A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative fre- 
quencies, x‘(y). (C) Changes in the average levels of cooperation, private trait expression, and 
pleiotropy, Z‘(y). (D) Differential fitness effects of loss-of-function mutations within the 
group. Parameters: s° = 0.95; K = 200; 4 = 0.0001; v = 0.01. The code required to generate this 
figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/ 
record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$4 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 

i= 4. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype g,, which actively expresses a private trait, z; = 1, and pleiotropy, z; = 1, but no 
cooperation, z; = 0. Growth of the group as its age, y, increases, is logistic, with a carrying 
capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, 
¢. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 
25% of the group. Note that the strength of pleiotropy has no effect on the within-group 
dynamics. (A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative fre- 
quencies, x‘(y). (C) Changes in the average levels of cooperation, private trait expression, and 
pleiotropy, Z‘(y). (D) Differential fitness effects of loss-of-function mutations within the 
group. Parameters: s° = 0.95; K = 200; u = 0.0001; v = 0.01. The code required to generate this 
figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/ 
record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$5 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 

i= 3. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype gu, which actively expresses a private trait, z° = 1, but no cooperation, z‘ = 0, or plei- 
otropy, z; = 0. Growth of the group as its age, y, increases, is logistic, with a carrying capacity 
K =200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, ¢. The 
vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 25% 
of the group. Note that the strength of pleiotropy has no effect on the within-group dynamics. 
(A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative frequencies, 
x“(y). (C) Changes in the average levels of cooperation, private trait expression, and pleiotropy, 
Z‘(y). (D) Differential fitness effects of loss-of-function mutations within the group. Parame- 
ters: s° = 0.95; K = 200; w = 0.0001; v = 0.01. The code required to generate this figure can be 
found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#. 
YjSBVurP2Uk. 

(DOCX) 





S6 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 
i= 2. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype g», which actively expresses pleiotropy, z; = 1, but no cooperation, zi, = 0, a no pri- 


vate trait, z° = 0. Growth of the group as its age, y, increases, is logistic, with a carrying 
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capacity K = 200. Dynamics are shown from left to right for 3 different strengths of pleiotropy, 
¢. The vertical dashed line in (A-C) represents the point at which mutant cell lineages make up 
25% of the group. Note that the strength of pleiotropy has no effect on the within-group 
dynamics. (A) Changes in genotype abundances, n‘(y). (B) Changes in genotype relative fre- 
quencies, x‘(y). (C) Changes in the average levels of cooperation, private trait expression, and 
pleiotropy, Z‘(y). (D) Differential fitness effects of loss-of-function mutations within the 
group. Parameters: s° = 0.95; K = 200; 4 = 0.0001; v = 0.01. The code required to generate this 





figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/ 
record/6367788#.YjSBVurP2Uk. 
(DOCX) 


$7 Fig. Within-group evolutionary dynamics of groups founded by cells with genotype 

i= 1. Within-group mutation selection dynamics are shown for a group founded by a cell with 
genotype g,, which does not actively expresses any of cooperation, z° = 0, a private trait, 

zi, = 0, or pleiotropy, z; = 0. Growth of the group as its age, y, increases, is logistic, with a car- 
rying capacity K = 200. Dynamics are shown from left to right for 3 different strengths of plei- 
otropy, ¢. The vertical dashed line in (A-C) represents the point at which mutant cell lineages 
make up 25% of the group. Note that the strength of pleiotropy has no effect on the within- 
group dynamics. (A) Changes in genotype abundances, n“(y). (B) Changes in genotype relative 
frequencies, x“(y). (C) Changes in the average levels of cooperation, private trait expression, 
and pleiotropy, Z‘(y). The differential fitness effects of loss-of-function mutations are not 
shown, because g; cells have no traits to lose. Parameters: s° = 0.95; K = 200; v = 0.0001; v = 
0.01. The code required to generate this figure can be found at https://github.com/euler-mab/ 
pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$8 Fig. Long-term evolutionary dynamics to illustrate a case where the evolution of cooper- 
ation depends entirely on pleiotropy. Dynamics are shown for the global population of cells 
over time, t, which encompasses many generations of cell groups. These dynamics, therefore, 
encompass both within-group and between-group selection dynamics. We show 3 strengths of 
pleiotropy, ¢, to capture 3 qualitatively different scenarios. Stronger pleiotropy is associated 
with more rapid and complete evolution of both pleiotropy and cooperation. (A) Changes in 
global genotype relative frequencies, x“(t). (B) Changes in the global average levels of coopera- 
tion, private trait expression, and pleiotropy, Z‘(t). In the absence of pleiotropy ¢ = 0, coopera- 
tion fails to evolve. (C) Average change in traits over a group lifetime, measured as the 
difference between the average trait values among groups aged y to those expected from their 
founding cell at birth. Parameters: s° = s* = 0.95; K = 500; = 0.0001; v = 0.01; A = 25. The code 
required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and 
https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$9 Fig. Private trait evolution does not favour pleiotropy and vice versa. We explored a 
model in which the cooperative trait was replaced by another private trait. Heatmaps show 
average trait values of 2 private traits and pleiotropy among the global population of cells 
(across all groups) at steady state in our revised model. Results are shown for 3 group sizes 
(increasing from top to bottom). Both private traits evolve to fixation under all parameter val- 
ues, but pleiotropy is never favoured. The dotted line marks the boundary between pleiotropy 
having no effect (control case) and pleiotropy having an effect on the outcome of mutations. 
Parameters: s° = 0.95; K = 200; u = 0.0001; v = 0.01. The code required to generate this figure 
can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/ 
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6367788#.YjSBVurP2Uk. 
(DOCX) 


$10 Fig. Within-group evolutionary dynamics of groups founded by 8 cells with all 8 geno- 
types represented. Within-group mutation selection dynamics are shown for a group founded 
by cells with all possible genotypes and phenotypes represented. Growth of the group as its 
age, y, increases, is logistic, with a carrying capacity K = 200. Dynamics are shown from left to 
right for 3 different strengths of pleiotropy, ¢. Note that the strength of pleiotropy has no effect 
on the within-group dynamics. (A) Changes in genotype abundances, n‘(y). (B) Changes in 
genotype relative frequencies, x“(y). (C) Changes in the average levels of cooperation, private 
trait expression, and pleiotropy, Z‘(y). (D) Differential fitness effects of loss-of-function muta- 
tions within the group. Parameters: s° = 0.95; K = 200; u = 0.0001; v = 0.01. The code required 
to generate this figure can be found at https://github.com/euler-mab/pleiotropy and https:// 
zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$11 Fig. Pleiotropy stabilises cooperation when groups are founded by 2 cells rather than 
one. We explored a model in which 2 cells are selected uniformly at random from the ancestor 
group to found each descendant group (rather than a single cell, which we assume in the main 
paper). This assumption lowers the expected relatedness in groups at the point at which they 
form. Heatmaps show average trait values among the global population of cells (across all 
groups) at steady state in our model. Results are shown for 3 group sizes (increasing from top 
to bottom). As when groups are founded by a single cell, pleiotropy is favoured when the 
strength of pleiotropy, ¢, is higher. The overall levels of cooperation are lower when groups are 
founded by 2 cells rather than one, but the evolution of pleiotropy still promotes the evolution 
of cooperation. The dotted line marks the boundary between pleiotropy having no effect (con- 
trol case) and pleiotropy having an effect on the outcome of mutations. Parameters: s° = s* = 
0.95; K = 200; = 0.0001; v = 0.01. The code required to generate this figure can be found at 
https://github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#. 
YjSBVurP2Uk. 

(DOCX) 


$12 Fig. Pleiotropy evolves to stabilise cooperation even when costly. We explored a model 
in which pleiotropy carries a cost on multicellular function. We assumed multicellular group 
function is reduced by a factor 1 — Zz), where ¢ is the cost of pleiotropy, and Z; is the average 
pleiotropy in a group. Heatmaps show average trait values among the global population of 
cells (across all groups) at steady state in our model. Results are shown for 3 levels of cost 
(increasing from top to bottom). Pleiotropy evolves even when it carries a 1%, 2%, or 5% cost 
to group functionality. When costly, it is especially favoured when groups are longer lived and 
therefore require mechanisms to limit the spread of noncooperative mutant lineages. The dot- 
ted line marks the boundary between pleiotropy having no effect (control case) and pleiotropy 
having an effect on the outcome of mutations. Parameters: s° = s* = 0.95; K = 200; u = 0.0001; v 
= 0.01; K = 200. The code required to generate this figure can be found at https://github.com/ 
euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 








$13 Fig. Pleiotropy evolves to stabilise cooperation even when group selection is weak. We 
explored a model in which the strength of group selection, s°, is varied. Heatmaps show aver- 
age trait values among the global population of cells (across all groups) at steady state in our 
model. Results are shown for 3 strengths of group selection (increasing from top to bottom). 
Stronger group selection favours both cooperation and pleiotropy, but strong pleiotropy can 
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help cooperation evolve even when group selection is weak, (s® = 0.05). The dotted line marks 
the boundary between pleiotropy having no effect (control case) and pleiotropy having an 
effect on the outcome of mutations. Parameters: s° = 0.95; K = 200; 4 = 0.0001; v = 0.01; 

K = 200. The code required to generate this figure can be found at https://github.com/euler- 
mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$14 Fig. The rate of gain-of-function mutations has little impact on the dynamics of pleiot- 
ropy evolution. We varied the relative rate of gain-of-function mutations to loss-of-function 
mutations, v. Heatmaps show average trait values among the global population of cells (across 
all groups) at steady state in our model. Results are shown for 3 gain-of-function rates (increas- 
ing from top to bottom). The gain-of-function ratio had a marginal effect on the evolution of 
pleiotropy, favouring slightly higher rates when v is higher. The dotted line marks the bound- 
ary between pleiotropy having no effect (control case) and pleiotropy having an effect on the 
outcome of mutations. Parameters: s° = s* = 0.95; K = 200; ¢ = 0.0001; K = 200. The code 
required to generate this figure can be found at https://github.com/euler-mab/pleiotropy and 
https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





S15 Fig. Strong pleiotropy can help rescue cooperation even in the face of high mutation 
rates. We varied the loss-of-function mutation rate, 4. Heatmaps show average trait values 
among the global population of cells (across all groups) at steady state in our model. Results 
are shown for 3 loss-of-function rates (increasing from top to bottom). Higher mutation rates 
disfavoured the evolution of cooperation, but pleiotropy still evolved at higher strengths of 
pleiotropy and for lower group sizes. The evolution of pleiotropy was associated with a stabili- 
sation of cooperation even at higher mutation rates. The dotted line marks the boundary 
between pleiotropy having no effect (control case) and pleiotropy having an effect on the out- 
come of mutations. Parameters: s° = s* = 0.95; K = 200; v = 0.01. The code required to generate 
this figure can be found at https://github.com/euler-mab/pleiotropy and https://zenodo.org/ 
record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$16 Fig. Lower mutation rates increases cooperation but pleiotropy still evolves. We varied 
the loss-of-function mutation rate, 1, using values 10-fold lower than the typical value used in 
the main text. Heatmaps show average trait values among the global population of cells (across 
all groups) at steady state in our model. Results are shown for 3 loss-of-function rates (increas- 
ing from top to bottom). With 10-fold lower mutation rates, we found that cooperation breaks 
down only for longer group life spans, 2. When cooperation is threatened by breakdown, how- 
ever, we still found that pleiotropy evolves to rescue it. The dotted line marks the boundary 
between pleiotropy having no effect (control case) and pleiotropy having an effect on the out- 
come of mutations. Parameters: s° = s* = 0.95; K = 200; v = 0.01; c = 0. The code required to 
generate this figure can be found at https://github.com/euler-mab/pleiotropy and https:// 
zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$17 Fig. Lower mutation rates increase cooperation but pleiotropy still evolves, even when 
pleiotropy is costly. We varied the gain-of-function mutation rate, 4, using values 10-fold 
lower than the typical value used in the main text, and assumed that the evolution of pleiotropy 
decreases group function by 2%. Heatmaps show average trait values among the global popula- 
tion of cells (across all groups) at steady state in our model. Results are shown for 3 loss-of- 
function rates (increasing from top to bottom). With 10-fold lower mutation rates and a cost, 
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we found that pleiotropy evolves most frequently when the mutation rate is higher and cooper- 
ation is under the greatest threat. The dotted line marks the boundary between pleiotropy hav- 
ing no effect (control case) and pleiotropy having an effect on the outcome of mutations. 
Parameters: s° = s* = 0.95; K = 200; v = 0.01; € = 0.02. The code used to produce this figure may 
be found at https://github.com/euler-mab/pleiotropy. 

(DOCX) 





$18 Fig. Reproduction at maturity hinders cooperation and has mixed effects on the evolu- 
tion of pleiotropy. We varied the age fraction of the expected life span at which groups reach 
reproductive maturity, a, and assumed that the evolution of pleiotropy decreases group func- 
tion by 2%. Heatmaps show average trait values among the global population of cells (across 
all groups) at steady state in our model. Results are shown for 3 reproductive maturity parame- 
ters (increasing the age at which maturity is reached from top to bottom). Increasing the age of 
reproductive maturity has a marginal positive effect on the evolution of pleiotropy for shorter 
life spans, but a marginal negative effect on cooperation. The dotted line marks the boundary 
between pleiotropy having no effect (control case) and pleiotropy having an effect on the out- 
come of mutations. Parameters: s“ = s* = 0.95; K = 200; v = 0.01; ¢ = 0.02. The code required to 
generate this figure can be found at https://github.com/euler-mab/pleiotropy and https:// 
zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 





$19 Fig. When mutation-driven breakdown of cooperation is a problem, pleiotropy stabi- 
lises cooperation. We modified the individual-based model of dos Santos and colleagues to 
compare what happens when pleiotropy cannot evolve (left) as compared to when it can evolve 
(right). We varied the length of the within-group growth phase k during which spontaneous 
mutants can arise and invade within groups, and the relatedness r at the point at which groups 
form. We follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evo- 
lutionary dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) 
Steady-state levels of all 3 traits under a when the length of the within-group growth phase is 
varied (x-axis) for r = 1. (C) Steady-state levels of all 3 traits when relatedness at the point at 
which groups form is varied for a within-group growth phase of k = 30. Other parameters: 
b=0.11,c=0.1, g=0.5, mutation rate u = 0.001, number of groups Ng = 1000. All plots are 
averages of 10 replicates. The code required to generate this figure can be found at https:// 
github.com/euler-mab/pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 
(DOCX) 


$20 Fig. Pleiotropy evolves to stabilise cooperation even when the loss-of-function muta- 
tion rate is 10x greater than the gain-of-function rate. We performed simulations using the 
mutation model used by dos Santos and colleagues in the production of their $5 Fig. This 
again reveals conditions where pleiotropy evolves and increases cooperation as it does. We var- 
ied the length of the within-group growth phase k during which spontaneous mutants can 
arise and invade within groups, and the relatedness r at the point at which groups form. We 
follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evolutionary 
dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) Steady-state 
levels of all 3 traits under a when the length of the within-group growth phase is varied (x-axis) 
for r= 1. (C) Steady-state levels of all 3 traits when relatedness at the point at which groups 
form is varied for a within-group growth phase of k = 30. Other parameters: b = 0.11, c= 0.1, 
g = 0.5, mutation rate v = 0.001, number of groups ng = 1000. All plots are averages of 10 repli- 
cates. The code required to generate this figure can be found at https://github.com/euler-mab/ 
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pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 
(DOCX) 


$21 Fig. Pleiotropy evolves to stabilise cooperation even when the loss-of-function muta- 
tion rate is 100x greater than the gain-of-function rate. We performed simulations using 
the mutation model used by dos Santos and colleagues in the production of their S6 Fig. This 
again reveals conditions where pleiotropy evolves and increases cooperation as it does. We var- 
ied the length of the within-group growth phase k during which spontaneous mutants can 
arise and invade within groups, and the relatedness r at the point at which groups form. We 
follow the evolution of a private trait, cooperative trait, and pleiotropy trait. (A) Evolutionary 
dynamics of all 3 traits for a within-group growth phase of k = 30 and r = 1. (B) Steady-state 
levels of all 3 traits under a when the length of the within-group growth phase is varied (x-axis) 
for r = 1. (C) Steady-state levels of all 3 traits when relatedness at the point at which groups 
form is varied for a within-group growth phase of k = 30. Other parameters: b = 0.11, c= 0.1, 

g = 0.5, mutation rate v = 0.001, number of groups ng = 1000. All plots are averages of 10 repli- 
cates. The code required to generate this figure can be found at https://github.com/euler-mab/ 
pleiotropy and https://zenodo.org/record/6367788#.YjSBVurP2Uk. 

(DOCX) 


Acknowledgments 


With thanks to Dave Queller and Arvid Agren for helpful comments on an earlier version of 
the main text, to Stu West and Miguel dos Santos for helping us to better understand their 
own model of pleiotropy and to Andrew Murray for a formative conversation on the evolution 
of pleiotropy. 


Author Contributions 

Conceptualization: Michael A. Bentley, Christian A. Yates, Gail M. Preston, Kevin R. Foster. 
Data curation: Michael A. Bentley. 

Formal analysis: Michael A. Bentley. 

Funding acquisition: Kevin R. Foster. 

Investigation: Michael A. Bentley. 

Software: Michael A. Bentley. 

Supervision: Christian A. Yates, Jotun Hein, Gail M. Preston, Kevin R. Foster. 

Writing - original draft: Michael A. Bentley, Kevin R. Foster. 


Writing - review & editing: Michael A. Bentley, Christian A. Yates, Jotun Hein, Gail M. 
Preston, Kevin R. Foster. 


References 


1. Nadell CD, Drescher K, Foster KR. Spatial structure, cooperation and competition in biofilms. Nat Rev 
Microbiol. 2016; 14:589-600. htips://doi.org/10.1038/nrmicro.2016.84 PMID: 27452230 


2. Fisher RM, Cornwallis CK, West SA. Group formation, relatedness, and the evolution of multicellularity. 
Curr Biol. 2013; 23:1120—5. https://doi.org/10.1016/j.cub.2013.05.004 PMID: 23746639 


3. Ghoul M, Griffin AS, West SA. Toward an evolutionary definition of cheating. Evolution. 2014. https:// 
doi.org/10.1111/evo.12266 PMID: 24131102 


4. Greig D, Travisano M. The Prisoner’s Dilemma and polymorphism in yeast SUC genes. Proc R Soc B 
Biol Sci. 2004; 271:25-6. https://doi.org/10.1098/rsbl.2003.0083 PMID: 15101409 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001626 June 3, 2022 31/34 


PLOS BIOLOGY 


Pleiotropic constraints promote the evolution of cooperation in cellular groups 





10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


Diggle SP, Griffin AS, Campbell GS, West S, a. Cooperation and conflict in quorum-sensing bacterial 
populations. Nature. 2007; 450:41 1-4. https://doi.org/10.1038/nature06279 PMID: 18004383 


Harrison F, Browning LE, Vos M, Buckling A. Cooperation and virulence in acute Pseudomonas aerugi- 
nosa infections. BMC Biol. 2006; 4:1—5. https://doi.org/10.1186/1741-7007-4-21 


Griffin AS, West SA, Buckling A. Cooperation and competition in pathogenic bacteria. Nature. 2004; 
430:1024-7. https://doi.org/10.1038/nature02744 PMID: 15329720 


Sandoz KM, Mitzimberg SM, Schuster M. Social cheating in Pseudomonas aeruginosa quorum sens- 
ing. Proc Natl Acad Sci U S A. 2007; 104:15876-81. https://doi.org/10.1073/pnas.0705653104 PMID: 
17898171 


Velicer GJ, Kroos L, Lenski RE. Developmental cheating in the social bacterium Myxococcus xanthus. 
Nature. 2000; 404:598-601. https://doi.org/10.1038/35007066 PMID: 10766241 


Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL. Solutions to the Public Goods Dilemma in 
Bacterial Biofilms. Curr Biol. 2014; 24:50-5. https://doi.org/10.1016/j.cub.2013.10.030 PMID: 
24332540 


Andersen SB, Marvig RL, Molin S, Johansen HK, Griffin AS. Long-term social dynamics drive loss of 
function in pathogenic bacteria. Proc Natl Acad Sci U S A. 2015; 112:10756-61. https://doi.org/10. 
1073/pnas. 1508324112 PMID: 26240352 


Raymond B, West §S, a, Griffin AS, Bonsall MB.. The dynamics of cooperative bacterial virulence in the 
field. Science. 2012; 337:85-8. https://doi.org/10.1126/science. 1218196 PMID: 22767928 


Cordero OX, Ventouras L-A, DeLong EF, Polz MF. Public good dynamics drive evolution of iron acquisi- 
tion strategies in natural bacterioplankton populations. Proc Natl Acad Sci U S A. 2012; 109:20059-64. 
https://doi.org/10.1073/pnas. 1213344109 PMID: 23169633 


Grosberg RK, Strathmann RR. The Evolution of Multicellularity: A Minor Major Transition? Annu Rev 
Ecol Evol Syst. 2007; 38:621—54. https://doi.org/10.1146/annurev.ecolsys.36.102403.114735 


Grosberg RK, Strathmann RR. One cell, two cell, red cell, blue cell: The persistence of a unicellular 
stage in multicellular life histories. Trends Ecol Evol. 1998; 13:112-6. https://doi.org/10.1016/S0169- 
5347(97)01313-X PMID: 21238226 


Hamilton WD. The genetical evolution of social behaviour. | & II J Theor Biol. 1964; 7:1—52. 


Foster KR. A defense of sociobiology. Cold Spring Harb Symp Quant Biol. 2009; 74:403-18. https://doi. 
org/10.1101/sqb.2009.74.041 PMID: 20375319 


Michod RE. Evolution of the individual. Am Nat. 1997; 150:S5-s21. https://doi.org/10.1086/286047 
PMID: 18811312 


Michod RE, Roze D. Cooperation and conflict in the evolution of multicellularity. Heredity. 2001:1—7. 
https://doi.org/10.1046/j. 1365-2540.2001 .00808.x PMID: 11298810 


Michod RE. Darwinian dynamics: evolutionary transitions in fitness and individuality. Princeton 
University Press; 2000. 





Dandekar AA, Chugani S, Greenberg EP. Bacterial quorum sensing and metabolic incentives to coop- 
erate. Science (80-). 2012; 338:264—6. https://doi.org/10.1126/science. 1227289 PMID: 23066081 


Wang M, Schaefer AL, Dandekar AA, Greenberg EP. Quorum sensing and policing of Pseudomonas 
aeruginosa social cheaters. Proc Natl Acad Sci U S A. 2015; 112:2187-91. https://doi.org/10.1073/ 
pnas.1500704112 PMID: 25646454 


Sathe S, Mathew A, Agnoli K, Eberl L, Kummerli R. Genetic architecture constrains exploitation of side- 
rophore cooperation in the bacterium Burkholderia cenocepacia. Evol Lett. 2019; 3:610-22. https://doi. 
org/10.1002/evI3.144 PMID: 31844554 


Foster KR, Shaulsky G, Strassmann JE, Queller DC, Thompson CRL. Pleiotropy as a mechanism to 
stabilize cooperation. Nature. 2004; 431:693-6. https://doi.org/10.1038/nature02894 PMID: 15470429 


Foster KR, Parkinson K, Thompson CRL. What can microbial genetics teach sociobiology? Trends 
Genet. 2007; 23:74-80. https://doi.org/10.1016/j.tig.2006. 12.003 PMID: 17207887 


Frénoy A, Taddei F, Misevic D. Genetic architecture promotes the evolution and maintenance of coop- 
eration. PLoS Comput Biol. 2013; 9:e1003339. htips://doi.org/10.1371/journal.pcbi. 1003339 PMID: 
24278000 


Chisholm RH, Connelly BD, Kerr B, Tanaka MM. The role of pleiotropy in the evolutionary maintenance 
of positive niche construction. Am Nat. 2018; 192:35—48. https://doi.org/10.1086/697471 PMID: 
29897798 


Sharp C, Foster KR. 2022. Host control and the evolution of cooperation in host microbiomes. Nature 
Communications, in press. 


dos Santos M, Ghoul M, West SA. Pleiotropy, cooperation, and the social evolution of genetic architec- 
ture. PLoS Biol. 2018; 16:1—25. https://doi.org/10.1371/journal.pbio.2006671 PMID: 30359363 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001626 June 3, 2022 32/34 


PLOS BIOLOGY 


Pleiotropic constraints promote the evolution of cooperation in cellular groups 





30. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


47. 


48. 


49. 


50. 


51. 


52. 


53. 


54. 


55. 


David CN, Campbell RD. Cell Cycle Kinetics and Development of Hydra Attenuata: |. Epithelial Cells J 
Cell Sci. 1972; 11:557-68. https://doi.org/10.1242/jcs.11.2.557 PMID: 5076361 


Foster KR. The sociobiology of molecular systems. Nat Rev Genet. 2011; 12:193—203. hitps://doi.org/ 
10.1038/nrg2903 PMID: 21301472 


Grafen A. Natural selection, kin selection and group selection. In: Krebs JR, Davies NB, editors. Beha- 
vioural ecology: An evolutionary approach. Oxford: Blackwell; 1984. p. 62-84. 


Paaby AB, Rockman MV. The many faces of pleiotropy. Trends Genet. 2013; 29:66-—73. https://doi.org/ 
10.1016/j.tig.2012.10.010 PMID: 23140989 


Stearns FW. One hundred years of pleiotropy: a retrospective. Genetics. 2010; 186:767—73. https://doi. 
org/10.1534/genetics.110.122549 PMID: 21062962 


Guillaume F, Otto SP. Gene functional trade-offs and the evolution of pleiotropy. Genetics. 2012; 
192:1389-409. https://doi.org/10.1534/genetics.112.143214 PMID: 22982578 


Classon M, Harlow E. The retinoblastoma tumour suppressor in development and cancer. Nat Rev Can- 
cer. 2002; 2:910-7. https://doi.org/10.1038/nrc950 PMID: 12459729 


Hickman ES, Moroni MC, Helin K. The role of p53 and pRB in apoptosis and cancer. Curr Opin Genet 
Dev. 2002; 12:60-6. https://doi.org/10.1016/s0959-437x(01)00265-9 PMID: 11790556 


Evan G, Lowe S, Cepero E. Intrinsic tumour suppression. Nature. 2004; 432:1-9. https://doi.org/10. 
1038/nature03098 PMID: 15549092 


Wagner GP, Zhang J. The pleiotropic structure of the genotype-phenotype map: The evolvability of 
complex organisms. Nat Rev Genet. 2011; 12:204—13. https://doi.org/10.1038/nrg2949 PMID: 
21331091 


Maynard-Smith J. Group selection and kin selection. Nature. 1964; 201:1145-7. 


Orr HA. Adaptation and the cost of complexity. Evolution (N Y). 2000; 54:13-20. https://doi.org/10. 
1111/j.0014-3820.2000.tb00002.x PMID: 10937178 


Milholland B, Dong X, Zhang L, Hao X, Suh Y, Vijg J. Differences between germline and somatic muta- 
tion rates in humans and mice. Nat Commun. 2017; 8:1-8. 


Werner B, Sottoriva A. Variation of mutational burden in healthy human tissues suggests non-random 
strand segregation and allows measuring somatic mutation rates. PLoS Comput Biol. 2018; 14: 
e1006233. https://doi.org/10.1371/journal.pcbi. 1006233 PMID: 29879111 


Tomczyk S, Fischer K, Austad S, Galliot B. Hydra, a powerful model for aging studies. Invertebr Reprod 
Dev. 2015; 59:11-6. https://doi.org/10.1080/07924259.2014.927805 PMID: 26120246 


Alivisatos AP, Chun M, Church GM, Greenspan RJ, Roukes ML, Yuste R. The Brain Activity Map Proj- 
ect and the Challenge of Functional Connectomics. Neuron. 2012; 74:970—4. htips://doi.org/10.1016/j. 
neuron.2012.06.006 PMID: 22726828 


Kim W, Racimo F, Schluter J, Levy SB, Foster KR. Importance of positioning for microbial evolution. 
Proc Natl Acad Sci U S A. 2014;111. https://doi.org/10.1073/pnas. 1323632111 


Bianconi E, Piovesan A, Facchin F, Beraudi A, Casadei R, Frabetti F, et al.. An estimation of the number 
of cells in the human body. Ann Hum Biol. 2013; 40:463-71. https://doi.org/10.3109/03014460.2013. 
807878 PMID: 23829164 


Mitri S, Foster KR. Pleiotropy and the low cost of individual traits promote cooperation. Evolution (N Y). 
2016; 70:488-94. https://doi.org/10.1111/evo.12851 PMID: 26748821 


Queller DC. Pleiotropy and synergistic cooperation. PLoS Biol. 2019:3—5. https://doi.org/10.1371/ 
journal.pbio.3000320 PMID: 31226105 


Korb J, Foster KR. Ecological competition favours cooperation in termite societies. Ecol Lett. 2010; 
13:754-60. https://doi.org/10.1111/.1461-0248.2010.01471.x PMID: 20584170 


Rankin DJ, L6pez-Sepulcre A, Foster KR, Kokko H. Species-level selection reduces selfishness 
through competitive exclusion. J Evol Biol. 2007; 20:1459-68. https://doi.org/10.1111/j.1420-9101. 
2007.01337.x PMID: 17584239 


Ostrowski EA. Enforcing Cooperation in the Social Amoebae. Curr Biol. 2019; 29:R474—84. https://doi. 
org/10.1016/j.cub.2019.04.022 PMID: 31163162 


Wielgoss S, Schneider D, Barrick JE, Tenaillon O, Cruveiller S, Chane-Woon-Ming B, et al.. Mutation 
rate inferred from synonymous substitutions in a long-term evolution experiment with escherichia coli. 
G3 Genes, Genomes Genet. 2011; 1:183-6. https://doi.org/10.1534/g3.111.000406 PMID: 22207905 


Crespi BJ. The evolution of social behavior in microorganisms. Trends Ecol Evol. 2001; 16:178-83. 
https://doi.org/10.1016/s0169-5347(01)02115-2 PMID: 11245940 


Harrison F, Buckling A. Hypermutability impedes cooperation in pathogenic bacteria. Curr Biol. 2005; 
15:1968-71. https://doi.org/10.1016/j.cub.2005.09.048 PMID: 16271876 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001626 June 3, 2022 33/34 


PLOS BIOLOGY 


Pleiotropic constraints promote the evolution of cooperation in cellular groups 





56. 


57. 


58. 


59. 


60. 


61. 


62. 


63. 


64. 


65. 


66. 


67. 


68. 


69. 


70. 


71. 


72. 
73. 


74. 
75. 


Azimi S, Roberts AEL, Peng S, Weitz JS, McNally A, Brown SP, et al.. Allelic polymorphism shapes 
community function in evolving Pseudomonas aeruginosa populations. ISME J. 2020; 14:1929-42. 
https://doi.org/10.1038/s41396-020-0652-0 PMID: 32341475 


Taylor TB, Johnson LJ, Jackson RW, Brockhurst MA, Dash PR. First steps in experimental cancer evo- 
lution. Evol Appl. 2013; 6:535—48. https://doi.org/10.1111/eva.12041 PMID: 23745144 


Lean C, Plutynski A. The evolution of failure: explaining cancer as an evolutionary process. Biol Philos. 
2016; 31:39-57. 


Nunney L. Lineage selection and the evolution of multistage carcinogenesis. Proc R Soc B Biol Sci. 
1999; 266:493-8. https://doi.org/10.1098/rspb.1999.0664 PMID: 10189713 


Nunney L. The real war on cancer: the evolutionary dynamics of cancer suppression. Evol Appl. 2013; 
6:11-9. https://doi.org/10.1111/eva.12018 PMID: 23396311 


Singh M, Boomsma JJ. Policing and punishment across the domains of social evolution. Oikos. 2015; 
124:971-82. 


DeGregori J. Adaptive oncogenesis: a new understanding of how cancer evolves inside us. Harvard 
University Press; 2018. 


DeGregori J. Challenging the axiom: does the occurrence of oncogenic mutations truly limit cancer 
development with age? Oncogene. 2013; 32:1869-75. https://doi.org/10.1038/onc.2012.281 PMID: 
22751134 


Serrano M, Lin AW, McCurrach ME, Beach D, Lowe SW. Oncogenic ras provokes premature cell 
senescence associated with accumulation of p53 and p16INK4a. Cell. 1997; 88:593-602. https://doi. 
org/10.1016/s0092-8674(00)81902-9 PMID: 9054499 


Schreiber RD, Old LU, Smyth MJ. Cancer immunoediting: integrating immunity’s roles in cancer sup- 
pression and promotion. Science (80-). 2011; 331:1565—70. https://doi.org/10.1126/science. 1203486 
PMID: 21436444 


Alpar L, Bergantifios C, Johnston LA. Spatially restricted regulation of Spatzle/Toll signaling during cell 
competition. Dev Cell. 2018; 46:706-19. https://doi.org/10.1016/j.devcel.2018.08.001 PMID: 30146479 


Russell JH, Ley TJ. Lymphocyte-mediated cytotoxicity. Annu Rev Immunol. 2002; 20:323-70. https:// 
doi.org/10.1146/annurev.immunol.20.100201.131730 PMID: 11861606 


Kajita M, Fujita Y. EDAC: Epithelial defence against cancer—cell competition between normal and 
transformed epithelial cells in mammals. J Biochem. 2015; 158:15-23. https://doi.org/10.1093/jb/ 
mvv050 PMID: 25991731 


C. QD, Eleonora P, Salvatore B, E. SJ. Single-Gene Greenbeard Effects in the Social Amoeba Dictyos- 
telium discoideum. Science (80-). 2003; 299: 105-106. https://doi.org/10.1126/science.1077742 PMID: 
12511650 


Schluter J, Schoech AP, Foster KR, Mitri S. The evolution of quorum sensing as a mechanism to infer 
kinship. PLoS Comput Biol. 2016; 12:e1004848. hitps://doi.org/10.1371/journal.pcbi.1004848 PMID: 
27120081 


Xavier JB, Kim W, Foster KR. A molecular mechanism that stabilizes cooperative secretions in Pseudo- 
monas aeruginosa. Mol Microbiol. 2011; 79:166-79. https://doi.org/10.1111/j.1365-2958.2010.07436.x 
PMID: 21166901 


Simon B, Fletcher JA, Doebeli M. Towards a General Theory of Group Selection; 2012. 


Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977; 81:2340- 
61. https://doi.org/10.1021/j100540a008 


LeVeque Ru. Finite volume methods for hyperbolic problems. Cambridge university press; 2002. 


Xavier JC, Patil KR, Rocha |. Metabolic models and gene essentiality data reveal essential and con- 
served metabolism in prokaryotes. PLoS Comput Biol. 2018; 14:e1006556. htips://doi.org/10.1371/ 
journal.pcbi. 1006556 PMID: 30444863 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3001626 June 3, 2022 34/34 


